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Chaptei 1 
INTRODUCTION 


When a solid object is placed in an otherwise uniform flow, vor- 
ticity develops in the vicinity of the solid boundary and in the wake 
as a result of the no-slip condition at the boundary. If the associ- 
ated Reynolds number is high enough, fluctuations in the flow appear 
as a result of the iiistability of the shear flow thus established. 
Fluctuations can develop in the wake as well as in the boundary layer 
and around the separation points and may lead to turbulence at rela- 
tively high Reynolds numbers. These flow fluctuations in turn, 
whether turbulent or not, induce a fluctuating pressure field on the 
surface of the body. Similarly, if the main stream Is turbulent or 
simgly unsteady, it gives rise, after being distorted by the presence 
of the body, to pressure fluctuations on its surface. The pressure 
fluctuations in question are of the so-called pseudo-sound type, and 
are due primarily to the instantaneous response of the pressure forces 
to fluctuations in the inertia and viscous forces in the neighborhood 
of the bjdy. 

The knowledge of the details of the pressure fluctuations occur- 
ring on the surface of a solid boundary interacting with a complicated 
flow field is of considerable importance in many practical problems. 

The sound radiated by solid surfaces in the presence of unsteady flows 
(Curie 1955; Ffowcs Williams and Hanwkings 1969) is partly due to a 
distribution of dipoles whose strengt’ depends upon the magnitude of 
the surface -pressure fluctuations. In many situations the contribution 
of these dipoles is the dominant one. In the recently recognized prob- 
lem of airframe noise which is due simply to the motion of the aircraft 
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body in the atmpsphere, the dominant part of the sound radiated results 
from the interaction of the unsteady flow structure in the vicinity of 
the aircraft with the aircraft body itself. Many other problems in the area 
of aerodynamic sound can be associated with similar types of interaction. 

In addition to the sound radiated away from a solid boundary supporting 
a turbulent flow, one often is concerned with die sound transmitted 
through that boundary, an aircraft fuselage for instance, in which case 
knowledge of the details of the pressure fluctuations acting on the sur- 
face is essential. Finally, in the area of structural design, whether 
applied to aeronautical problems or more generally to problems involv- 
ing the interaction of a structure with a fluid flow, it is very impor- 
tant to know the intensity and the scales of the surface-pressure fluc- 
tuations in order to determine the loads on the structure. 

Practical problems of the type described abo^^e r^re usually very 
complicated and not amenable to direct analytical treatment. However, 
as is common in similar si.tiations, one begins by decomposing the given 
problem into several simpler ones, in the hope that after these simpler 
problems have been solved, their combination in one way or another 
would reveal in an approximate manner the general behavior of the solu- 
tion of the initial problem. Unfortunately in the type of situations 
we are concerned about, even after such a decomposition is made, and 
if the interesting features are to be preserved, one is still faced with 
complex phenomena such as separation and transition to turbulence, and 
to determine the details of the pressure fluctuations of interest via 
a complete solution of the problem is beyond the power of present day 
analytical techniques. In such circumstances one is compelled to seek 
a physical understanding of the phenomena involved in terms of cause 
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and effect in such a way as to identify each interacting component and 
to be able to fomulate its role mathematically. It is very important 
to note however that such distinction between cause and effect is hardly 
unique when nonlinear processes are involved. For instance, in the case 
of interest here and which involves a solid body interacting with a com- 
plicated flow field, perhaps the only unambiguous way of describing the 
physics of t'vi i.r^rblem is by writing the general Navier-Stokes equations 
and the associated boundary conditions. The equations themselves are 
simply the expression of the balance of the inertia, perssure, and vis- 
cous forces while the boundary conditions describe the presence of the 
body and the nature of the incoming flow. In such a case the value of 
any causal theory would depend upon its usefulness in rendering the 
problem tractable mathematically. The Acoustic Theory of Airflows, 
now forgotten, of Shaw (1949) is one example where the theory failed 
to lead to an analytically tractable formulation. In that theory the 
vortex shedding behind an airfoil in a uniform flow was regarded as con- 
trolled by pressure fluctuations on the surface of the airfoil. If, on 
the other hand, we regard the pressure fluctuations occurring on the 
surface of the body as the result of the interaction of an already 
existing and known velocity fi.eld with the geometry of the body, the 
Navier-Stokes equations lead under certain conditions to a linear equa- 
tion for the pressure of the Poisson type whose forcing term depends 
solely on the velocity field. The formal solution of such an equation 
is in the form of convolution product involving the forcing term ard an 
auxiliary function which depends exclusively on the geometry of the 
body, depicting in this manner the separate roles of the different ele- 
ments in the problem. Since it is rather the pressure field that we 
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are seeking here and since, in most cases, the velocity field is muif 
accessible to measurement or experimental techniques, the above ap]iroach 
is the most promising one, and consequently will be followed in the pre- 
sent work# 

Studies on the pressure fluctuations within a turbulent flow using 
the Poisson equation mentioned above began as early as 1948 with the 
work of Heisenberg, Obukhov, and Batchelor (1951). These early studies 
were concerned with the particular case of an unlAOi^ted field of homo- 
geneous and isotropic turbulence and used extensively the results of 
the theory of Homogeneous Turbulence which was developing rapidly at 
that time, Kraichnan (1956a), in prepa ation of a subsequent work on 
the case of a turbulent boundary layer introduced a scale of anisotropy 
into the problem, and studied the effect of a mean shear by just super- 
posing that shear on the initial field of homogeneous turbulence. In 
his second paper (Kraichnan, 1956b) he took into account the presence 
of an infinite plane boundary in an attempt to simulate the problem of 
the turbulent boundary layer. Although his work laid the foundation 
for subsequent works on the subject, his analysis was characterized by 
many artificial models which were only kinematically possible. Lilley 
and Hodgson (1960), and Lilley (1963) followed up on the subject using 
a different approach which took into account some of the physical pro- 
perties of the flow known from velocity measurements in the turbulent 
boundary layer, but Inevitably encumbered the analysis with many approx- 
imations and crude assumptions. Nevertheless, one interesting result 
emerged from their work which suggests that the wall pressure fluctua- 
tions are mainly due to contributions from sources located in the outer 
mixing region of the boundary layer where large eddies are swept in the 
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flow direction with relatively high speeds. Corcos (1964) elaborated 
further on this point in an interesting paper devoted to a discussion 
of the measurements of the statistical properties of the wall pressure 
and the associated difficulties. 

The theoretical work on the wall pressure fluctuations beneath a 
turbulent boundary layer, some of which was briefly described above, 
was motivated b) the great practical interest in that type of flows 
and the numerous experimental studies on the subject which appeared 
during the same period of time. But, due to the inherent difficulties 
of the structure of wall turbulence, the line of analytical approach 
based on the Poisson equation for the pressure was somewhat frustrated 
and as a result was abandoned for several years. Recently, however, 

Panton and Linebarger (1974) using that same equation attacked the prob- 

4 ^ 

lem with the help of the computer and obtained some interesting new re- 
sults. In their work the flow direction wavenumber spectrum of the wall 
pressure was expressed as a five -fold multiple integral which was inte- 
grated numerically using the Monte Carlo method, in throe different regions 
of the boundary layer, in order to assess the contribution of each region. 

Except for the case of a turbulent boundary on an infinite plane 
surface, the pressure fluctuations occurring on the surface of bodies 
of arbitrary shape in the presence of an unsteady flow have not at- 
fractecl the attention of theoretical workers. On the other hand, exam- 
ples where large-scale flow fluctuations not of the boundary layer type 


One exception is the work by Hunt (1973a) who considered the interaction 
of a bluff body with an incoming turbulent flow (Hunt, 1973b) on the 
basis of the theory of 'rapid distortion of turbulence' developed by 
Batchelor and Proudman. 
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occur In the vicinity of compact solid bodies, and CimsequeiUly indiic 
pressure fluctuations on their surfaces, are abundant and of considerable 
practical interest. The case of a cylindrical body lnter«acting with its 
own wake may be the prominent example and constitutes the sub ject of the 
present work. Some aspects of this research have been presented previously 
by the authors at the Second and Third Interagency Symposiums on University 
Research in Tran^sportatlon Noise (Karamcheti and Ayoub 1974; Ayoub and 
Karamcheti 1975'i . 

In the next chapter the problem of determining the characteri s tiCvS 

of the pressure fluctuations, induced on the surface ot a cylinder by the 

fluctuating wake hehind it is formulated. The flow in the wake is assumed 

incompressible and homogeneous in a direction along the span; the condition 

« * 

of incompressibility is discussed at die end of the chapter. A fonnal 
solution is then derived which relates the unsteady surface -pressure field 
to the velocity field In the w»ake. This allows certain general results to 
be derived, which are independent of the cylinder shape or the value of 
the Reynolds number. Although the next chapter is devoted primarily to 
the mathematical aspect of the problem, vsome divscussion of the physical 
aspect of the flow in the wake Is included, either to guide the mathematics 
or to interpret some new results. When this is done it Is generally in 
connection with bluff cylinders. 

In Chapter 3, the case of the circular cylinder is examined in de- 
tail, in light of the theory presented in Chapter 2. 
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PKOBLKM FORMIUATION ANI) GFI^.RAL ANALYSIS 

3 ♦ 1 S tatement ol the Prut I 

We consider an iufiu*vely long 
cylinder of constant but arbitrary 
cross-section placed in a wind which, 
far upstream, is uniform of magnitude 

V and in a direction normal to the 

■10 

cylinder axis. A system of orthogonal 
coordinates (OXj , Ox.^ , nx.j) as shown 
in the ) iguro, is set up. Let JU and 
D designate respectively the viscositj^ 
and de. 4 ty of rlie fluid which is as- 
sumed incompressible, and d a char- 
acteristic dimension of the cross-section. We assume that the value of 
the Reynolds number, which is defined by Re = (0 U^d)/u , is such that the 
perturbed flow is unsteady and three-dimensional, and also that the flow 
is homogeneous in the axial direction and of vanishing mean velocity compo- 
nent in that direction. The homogeneity in the axial direction implies 
that all mean values evaluated at a point in space are independent of the 
x^-coordinate , while cross-correlations at two points (x^ , , x^) and 

, x‘ , x') are functions of (x^-x*) and not of x« and x' sepa- 

i 2 ^ j j j J 

rately. We further assume that all fluctuating quantities are stationary 
functions of time, and therefore all mean values will be taken with respect 
to time according to the formula: 
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-!^r 

7(t7 lim ^ i7t) dt 

T ♦ CO U 

_T 

where the bar over a quantity is used to designate the mean of that quan- 
tity. It is proposed to study the eharncteristics of the pressure fluetua 
lions induced on the surface of the cylinder by the flow unsteadiness. 


1,1 Governing Kquations 

The flow being incompressible the governing Navier-Stokes ei|uations 
can be written in the fora: 



♦ 

where , v.^ , nre the three components of the velocity vector v , p 

the ]>ressure, t the time, u the kinematic viscosity defined by p 

and where the summation convection over repeated subscripts has been used. 

Equations (2.1) are to be vsupplemented by the no-slip condition at the 

solid boundary and appropriate conditions on the velocity and the pressure 

at infinity. By applying the operator on the first equation of (2.1) 

i 

and using the second, we obtain the following equation for the pressure: 


n 


-2 o“(v.v.) 

" P = _p LJ_ 

ftx . 9x . 5x . Bx . 


( 2 . 2 ) 


11 1 j 

If we now express the flow field as composed of two parts, one a steady 
mean part and the other a time dependent fluctuating part in the following 
manner: 
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Vj(x, " V. (Jt) + vl^(x, L) ; i = 

p(x,t) = p(x) + p' (x.t) 

wo obtain by .substituting the.se expressions Into (2.2), taking the tiiiic 
average, and .sub.straetlng the resulting equation from (2.2) itself, the 
f o 1 1 owing equa 1 1 on : 



wliere 


and 



q(x, t) 





1 .1 


(2.4) 


(:>/>) 


q. .(x,t) ” p(2v.vl -viv' ) . (2,6) 

1 .1 I I J I 

z\gai.n, 11* we su1)«t:irule (2.3) into the I’i.rst equation of (2,1), t.rdu^ 
the time average, and saibs tract the resulting equation from the original 
one, we obtain: 


P 





-'vl dv! o‘>I 

y I i. ^ y [ q. j/ 

1 ox . I ox . ox , ox . 
.) ' I .1 1 


(2.7) 


On the surface S of the cylinder assumed rigid, the mean and fluctuating 
arts of the velocity vector vanish separately and we have using (2.7): 


n. 

1 ox. 
1 


U n. ~ 

1 . ox . 5x . 
.1 J 


CMO 


where n is the unit vector (n^ , n^ , u) along the inward normal to S . 
For large values of |x| in a plane = constant , the mean velocity ap- 
proaches the constant value while the velocity fluctuation components 
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1 


decrease to zero like l/(x,“ + x,j ) inside the wake^ and like 
2 2 

l/(Xj^ +x., ) outside it. The behavior outside the wake can be seen as 
follows. The flow far upstream being considered uniform, the unsteady part 
of the vorticity vector is generated in the neighborhood of the cylinder, 
and decays at large distances from it. The irrotational velocity fluctua- 
tions outside the wake can then be thought of as induced by a concentration 
of vorticity in the near-wake region. But we know that the velocity in- 
duced by a distribution of vorticity (in an incompressible flow), which 
decays satisfactorily in all directions in space, behaves like l/jx| 
at large distances (see Batchelor 1967 § 2.9), Therefore since the wake 
is assumed to extend to infinity in the + x^-direction and is s taLii.tically 
homogeneous in that direction, the square >1;^ vc-iucit? fluctuation- m [ - 
side the wake behaves like. 


4-co 

r 



from which the behavior stated above follows. Referring to (2.7) we can 

, 2 

see that outside the wake decays at least like 1/ (x^ 

2 2 

and p^ at least like inside the wake and at large 

distances downstream the flow has reached the equilibrium stage and the 

pressure can be assumed to have approximately the value it has out- 

2 2 3/4 

side, i.e., p* decays like 1/ (x^ "^^2^ 


^Wsteady or turbulent wakes behind infinitely long cylinders remain so at 
arbitrarily large distances downstream. This is due to the fact that the 
Reynolds number based on the width W of the wake and a typical turbulent 
velocity is independent of the downstream coordinate 

and (see Landau and Lifshitz 1959 page 139). 
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Tho flow boinf? coitsidered homogeneous in the axial direction we per- 
form a generalized Fourier decomposition in tha* direction according to 
the relations: 


p' (x,t) 


4- on 4-00 

/ P(X,k,t)e '’dk , P(X,k,t) - / p'(x,t)e dx^ , 


v| (x,t) 


hf 


V^(X,k,tle ’dk , V^(X,k,t) 




i = 1,2,3 

+ ikx “ -ikx 

q(x,t)=j- / Q(X,k,t)c ^dk , Q(X,k,t) = r q(x,t)e dx. 


where X=(X^,X,j) is the position vector in the plane normal to the cy- 
linder, coinciding with the normal projection of x on that plane, k 
is a real variable, representing the wavenumber in the -ection, and 

P , Vi , Q are the complex Fourier amplitudes associated with k of the. 
physical quantities p' ,v!^,q . In the transformed vaiables equation 
(2.4) becomes: 


5"P , 2 

" 1’ - -Q • 

1 1 


( 2 . 10 ) 


Since n= (nj,n2,0) , Eq. (2.8) can be written explicitly as: 


, Bp' . Bp's 1 . 2 . 

1 2 j j J J s 


or since v| , v^ and consequently the derivative of vj and v^ 

vanish on S , as: 
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(n ^+n 

Bx, ^"2 Bx/ 


uCnj^C 


0 V, 


3“v' 


^*"ir * * *”V * 


1 ^ 0 "'^'9 

+ — ) +n ( =•+ ^)] 

Bx^ Bx2“^ Bxj“ Bx„“ 


Upon taking the Fourier transform with respect to x^ » and noting that 
n^^ and n^ are independent of x^ » we obtain; 


(n 


Bp 


1 BX 


+ n 


BP 


1 


2 BX 


) = U[n, ( 


-2 .2 

Vi B V, 


2 's 


1 '• 9 

'^2 


n ) C Q 

“ BX, 


9 ? 

^‘■V, B^V, i 


^+- 


4)^' 


BX. 


\s 


which can be written in a concise form as: 

9 


n 


BP 
1 BX. 




- V , I 

"" ^’^i BX.oX. ^ I 

J 1 ( S 


( 2.111 


The condition on P at infinity is the same as that on p' , i^e,, both 

1 


the real and imaginary parts of P vanish like 


X 


3/2 


for sufficiently 


large values of |xi . Consequently, P(X,k,t) appears to be the solution 
of the two-dimensional boundary value problem defined by equation (2.10), 
the boundai'" couvlition (2,11) and the condition at infinity stated above. 
This prob’’ is reproduced below for convenience: 


2 2 

V P-k P = -Q 


n.VP p = Un.VV 

L» 


( 2 . 12 ) 


.X 


3/2 


for X “♦ OD 


where uhe operators V and 7 are respectively the gradient and Laplace 

operators in two dimensions, C the line intersection of S with the 
plane (X^,X^) , and V the vector ' 

The differential equation in (2.12) is that which occurs in problems 
of steady diffusion with absorption proportional to the c..; ,.entration. 
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Henco for a given value of k , P con be interpreted os the steady 
temperature distribution in a two-dimensional domain of thermal conduc- 
tivity equal to 1, bounded internally be C and containing a source dis- 
tribution of heot, of density 4Q (amount of heat generated per unit area 
per unit time). Across C heat flows at a prescribed rate and the medium 
absorbs heat in proportion to the local temperature at the constant posi- 

O 

tive rate k"' . P is therefore the solution of a well -posed boundary 
value problem and hence is uniquely determined, if the boundary C 
satisfies certain conditions of regularity, by the data Q and the bound- 
ary values jLin.V-V . Note that the boundary condition on n.VP is not 
arbitrary in our problem but is compatible with the differential equation 
since they are both derived from the same Navier-Stokes equations. Hox^- 
ever, for the purpose of determining the pressure from a knowledge of the 

-4 

velocity field the value of n.VP on C can be considered as prescribed. 
The interpretation of P in terms of a temperature distribution is useful 
in assessing the relative importance of the different spectral componerts 

-4 2 “'* I 

of the pressure once some details on the distribution of Q and n.V V ^ 
for different values of k are known. 

2.3 Formal Solution . 

A formal integral representation of the solution of (2.12) can be 
derived using a Greenes function G(x|Y;k) defined by the following 

9 9 

V“G - k G = -6 (X-Y) 

"•®lc ■ 0 

G 0 as |x| 00 

"The fact that the coefficient of absorption is positive insures the 
uniqueness of the solution. 


) (2.13) 

J 
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where the differentiations are with 
respect to X , and Y^CY^jY^) . It 
can be ear.ily seen from the physical 
interpretation given in the preced- 
ing parapraph that (2, 13) admits a 
unique solution. Multiplying the 
first equation of (2.12) by G and 
the first equation of (2.13) by P , 
subtracting the resulting two equa- 
tions and integrating over the area 
bounded by C and a large cir- 
cle C of radius R centered at 
K 

the origin we obtain: 




(G7“P - P7*^G)dX = 


"R 



p6(X-Y)dX 
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where dX is an element of area. Using the Green's formula and Integra - 
ting the term containing O(X-Y) we obtain: 


P(Y,k,t) = 


I 

"r 


G QdX + 


/ 


(Gn.7P - Pn.7G)ds(X) 


(2.14) 


C+C 


R 


where n is the exterior normal to C and , and ds(X) is a line ele- 
ment of C and . It is easily seen from the behavior of G and P 

at large distances that as R goes to infinity the integral over 


vanishes. Furthermore, since n.7G =0 , and n.7P =/in.7 V ^ , re- 

f G j G I L 

lation (2.14) becomes: 
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P(Y,k,t) 


= Jf G(x1Y; 


k) Q(X,k,t)dX + 


U yG(xlY;k)n.V^V(X,k, 


t)ds(X) , CMi) 


where E is the region exterior to C extending to infinity. 

Equation (11.15) expresses the fact that at a given time t and at 
any point in space or on the cylinder having the coordinates and Y.^ 

in the Xj - and x^-directions, the component of the pressure associated 
with the wavenumber k (in the axial direction ) is determined uniquely 
by the corresponding Fourier components of the velocity fluctuations as 
distributed throughout the unsteady region in the (X^,X^) plane, and also 
by those pairs of components whose wavenumbers sums are equal to k . The 
contribution of the latter is due to the fact that q contains terms 
v.’hich are quadratic in the velocity fluctuations, and as can be easily 
seen the k-Fourier -component of the product instance, is no- 

thing but: 

4 * 03 


2n f V. (X,k,t)V.(X,k-k',t)dk' 


The surface integral in (2.15) can be thought of as the contribution to 
the pressure fluctuations by the Inertia forces and the integral over C 
as that of the viscous forces." Some details of the structure of the 
forcing term Q will be given in a later section. In the following we 
turn to the discussion of some practical considerations concerning the 


In fact the effect of the viscous forces is also included in the integral 
over X because these forces are responsible for maintaining a particular 
form of the inertia forces. However, the above distinction is made only 
in consistence with the approach adopted and explained in the introduction. 
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Greenes function G and the detemination of the total forces acting on 
the cylinder, 

2,4 The Greenes Function - Total Lift and Drag for the Case of a Circular 

Cylinder , 

Due to the linear character of the Poisson equation which governs the 
pressure , we were able to perform a spectral decomposition in the 
axial direction, so that each spectral component of the pre sure can be 
studied separately as the solution of a well-posed boundary value problem 
in two dimensions. That such a decomposition is meaningful physically fol- 
lows from the fact that the flow in general is homogeneous in the spanwise 
direction and is composed of two distinct parts: one is discrete and homo- 
geneous (or periodic for certain values of the Reynolds number) in that 
direction; the other is random in nature and can be assumed homogeneous in 
the same direction, in view of the two-dimensionality of the obstacle and 
of the flow upstream, A detailed discussion of this and other features of 
the flow behind a typical two-dimensional body (the circular cylinder) will 
be given in the next chapter. However, in order to guide our discussion 
of the Green's function in this section, a few words on the origin and the 
relative importance of these two components of the flow arj appropriate. 

When a bluff cylinder is placed in a uniform flow at relatively high 
Reynolds numbers, two boundary layers form on the sides of the cylinder, 
separate from its surface and tend to interact with each other to form the 
wake behind it. In this process two different phenomena take place inde- 
pendently. The first is the transition in the individual separated shear 
layers, which leads to small scale fluctuations and eventually to turbu- 
lence. The second is the direct interaction of these two shear layers. 
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carrying vorticity of opposite wsigns, which results in the formation of 
the large scale vortices characteristic of wakes behind blufl bodies at 
almost all Reynolds numbers. The random fluctuations have scales of the 
order of the thickness of shear layers emanating from the cylinder sur 
face and contain only a fraction of the wake energy, while the discrete 
component has length scales in the axial direction ranging from 2 to 30 
times a typical dimension of the cross-sect ion and carries most of the 
energy in the wake. These two compomuits which result from two indepen- 
dent phenomena in generally two different ranges of scales and frequencies 
preserve their separate identities for distance downstream and do not 

interact with each other in the near-wake region (see for instance Roshko, 
1954). With this in mind, it ivS clear that by studying the distribution 
of P on C for a given value of k , we evaluate the contribution to 
the pressure fluctuations on the cylinder, of a physically identifiable 
component of the flow. The usefulness of the interpretation of V in 
terms of a temperature distribution is likewise easily seen from the above 

discussion, for the wavenumber k is of the order of an inverse length 

2 

scale in the spanwise direction, and k is the rate of absorption in the 
analogue problem. 

In the form (2.15) the solution of the problem requires, among other 
things, a knowledge of the Green*s function G which satisfies (2.13). 

The usual method of finding G consists of writing it as follows: 

G(xlY;k) == G^CxlYjk) +G^(xlY;k) , ( 2 . 16 ) 

Here we are implicitly assuming that the contribution from the quadratic 
terms in the forcing function can be neglected. We will do the same in 
some of the discussions of §2.6 and §2.7. The role of these quadratic 
terms is discussed in Appendix A. 
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where Gy(Xly;k) is the free-space Green's function satisfying the first 
equation of (2.13) and the boundary condition at infinity, and Gj is 
a .solution of the corre.s ponding homogeneous equation, whidi is regular 
everywhere outside C and such that the sum G^ +Gj^ satisfie.s the second 
boundary condition in (2.13). G^^ is given by (see Stakgold, 1968): 

GQ(xlY;k) = •^Kjjrlk(X-Y)l] , (2.17) 

where i.s the zero-order modified Bessel function of the second kind. 

This function is tabulated (Abramowitz and Stegun, 1964), and has the fol- 
lowing asymptotic behavior for large and small values of the argument: 


K„(a) ~ 

/tT" -Oi 
' v -2a ^ 

for 

KjjCa) ~ 

' -logo: 

for a 0 


As for Gj^ , it satisfies the following problem: 

V^Gj^ - k^Gj^ = 0 

H.VGJ = -~n.7KQrlk(X-Y)l] 

1 c c 

0 as lx| 00 






(2.18) 


(2.19) 


For certain particular shapes of the boundary C , (2.19) can be solved 
by the method of separation of variables. The solution obtained is usu- 
ally in the fom of an infinite series of cylinder functions. When C is 
a circle centered at the origin, separatioii in the polar coordinates r 
and 0 yields an infinite sum of terms involving the circular functions 
e^^^ and the modified Bessel functions of the second kind K (|l^j) , 
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whosi- ;is\Tnptotic bc'havlor Tor large values of the argument is similar lo 
that of Kjj . Approximate solutions can then be found in some circum- 
stances depending on the value of k . Another way of approaching the 
problem is by using the free-spacc Green's function alone in (2.14), and 
instead of (2.15) the resulting relation becomes: 


P(Y,k,t) = G|j(X! Y;k) Q(X,k,t)dX+At J" G^j(xlY;k)n. V"V(X,k,t)ds(X) 

r 


- J P(X,k,t)n.7Gjj(X|Y;k)ds(X) 


( 2 . 20 ) 


2 "^ 

where is given by (2,17) • If Q and n.7"V are sssumed known this 

t'quation is an iiihoTnogeneous integral equation of the second kind in the 
curvilinear domain C , to v’hich in principle several approximate and nu- 
merical techniques can be applied to yield a solution for P on the sur- 
face of the cylinder. 

The methods of solution comtempla ted above are all approximate and 
have as a goal the determination of P , for a given value of k , as a 
function of position on the boundary C . Although it is usually desir- 
able to have such a solution it is often sufficient in practice, to ob- 
tain expressions for certain parameters, related in a global way to the 
detailed solution P . Among such parameters perhaps the most important 
ones, in the case of a bluff cylinder, are the total lift and drag forces 
acting on the cylinder. Using (2.20) we derive below such expressions for 
the case of a circulor cross-section. 


If such a course is to be followed (2.20) requires some modification un- 
less the dependence on time is simple, nevertheless in its present form 
Eq. (2.20) is useful as we will see below. 
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Lot r ,0 designate the polar coordinates of X and c , ; those 

d 

of Y , d/2 being the radius of C . We multiply (2.20) by -sin ;-d'. , 

after setting p = d/2 , and integrate from to r=2Ti , to obtain: 

2tt 

L(k,t) y~ -P(Y,k,t) sincp|dfp 
0 


2tt 

= J Gjj(xi\\k) Q(X,k,t)dX} 

0 r: 


2TT 2";T 

-i^) U y" sin'-pCj/" G^j(x!Y;k)n.V^V(X,k,t)d0 } d- 


0 

2tt 


0 

2rr 


'f sln,[ f P(X,k,t)n.7G^j(xlY;k)d0 ] df[),(2.2 
0 0 


1 ) 


where L(lcjt:) designates the fluctuating lift per unit span associated 
with the wavenumber k* , and evaluated at time t , We now set: 


2rr 

Qs(k,t) = GgQdxldcp 

0 X 


Sg(k,t) 



2tt 



G^n.y^^VdS ] dcp 




> ( 2 . 22 ) 
J 


and noticing that n. VG (xj Y;k) = - -r-^ (xlY;k) j 

^ |r=d/2 

changing the integration in the last term of (2.21): 


we 


obtain, 


after inter- 


The spanwise correlation length of the fluctuating surface-pressure at a 
particular frequency is found experimentally to be a function of the azi- 
muthal angle cp in general. The spectral decomposition in the x^-direc- 
tion has therefore the advantage of removing the ambiguity associated with 
the choice of a characteristic spanwise length scale of the fluctuating 
lift (or drag) at a given frequency. 
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J,ik, t ^ 


-f S ,{k,tU(“f /" P(X,k,t){^y^ G„(X|Y;k)sin 
0 0 


Tlu- integriil insidi* tlu* brackets In the 

last term of which we denote by 

1 , is evaluated at constant r and for 
* 

Y i*n the circle C , Now since 

Y;k) is a function t>l |X-Yj only, 

when Y Is on the circle, it becomes a 

i unction G t*r-0,r;k) of r and (r;i-0) 

0 

only. Therefore the integral in question 
can be written as: 
ll’'T 


1 


j' G||{';i-0,r;k) sin^ndto 


■d d0 
r d/ ' 

t;’..' 1 



or, after making the change of variable 'v='-p-0 , as: 

;'-T-0 


1 = 


G (v,r ;k) sin(v + 0)dv 


-0 

0 


2tt 


(y> t ;k) sin(v + 0)dv + J' (v,r ;k) sin(v + 0)dv 


-0 


2tt-0 


Gq (V , r ; k) s in (“V + 0) dv 


2rr 


Since G^(v + 2TT,r;k) = Gjj(V,r;k) the first and third terms in the above 
expression can be easily seen to cancel, and we have: 
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= co.s S J' G,|l,v,r;k) sinydv + siii 0 J GjjCy.r ;k ) i ds v dv 


lU’nce, C-’.-'k) bfccxnt'n: 


HkiL) = Q^(,k, t )+S^ (k , t 1 “(4) J (J,,^'V,r;k) sin V dvl ^ J l'(X,k , t i cv; 0 d0 

n n 

2tt ;>tt 

+ Ct" / C’',,lV,r';kjcos vdy} / P(X,k, t )sin 0d0 ] 

J “ J 


11' we define D(k,t) to be the I'liietuatinR droR )ier unit spim associated 
with the wavenumber k , we have: 


l)(k,t) -/ P(X,k,L)cos 0^d9 ; X on C 


and final Iv: 


L(k,L) = Q^(k,t) +S,,(k,t) ~ / (yy,r;k)slnvdy } jl)(k,t) 


+ / G.,(y,r ;k)cos y dy } ,L(.k,t) . (2.P.4) 


Similarly if we multiply (2.20) by -cos T — d^"* , and integrate from cr=0 
to rp = 2tt we obtain following the same procedure: 


l)(k,t) = Q_(k, t) +S (k, t) + — / G (y,r;k)cos y dv ] ,D(k,t) 


2" ;'r J “O' 
0 


f G (v,r;k)sinydy] ,L(k,L) 


( 2 . 2 ')'» 
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wht‘r»* Qj, ami S^, an- dfl incd as in (2,22) by replacing sin V by 
cos . Kqs. and (1*. 2'j) art? two independent algebraic equations in 

I. and n , and hence yield the following expressions for the lift and 
drag forces I-(k,t) and DCk,!) ; 


L(k,L) 


fq,,(k,ti +Sj,(k,t)}{l -G,j|.(k)} +{Q^(k,t) +S^(k,t)' ^jj.(k)^ 


{1-G^^(k)f -t-G.^yCk)' 




os 


{ Q ( k , t ) f S, . (k , t ) } { 1 -G (k; } -C Q„ (k , t ) + S t'k , t) ] G,,_ (k) 
^ 

1 2 


{l-G'j^(k)] +Gf)j.(k)' 


where 


I’rr 


n|,,.lk) = 7 {^ f G (v,r;k) cosvdv] 


J. ' -?r (; 

0 


r- d/; 


*!»("» • ftv J »«-yv5 


An import ;ml: foaturo of the formal solution and the OMprossion' 

C'.lMvi and (?, ,27) is that the Integrals over Involved converge very 
rapidly. These integrals which in principle should be integrated over 
the entire rt^gion wht^re Q is different fr<Mu zero, can be integr ;t:ed in 
a practical problem over a limited region around the point Y , due to 
the asymptotic behavior of G and for large valuers of [x-Yj . If 

P is to be evaluated on C the domain of integration is generally con- 
fined to the near-wake and is smaller the larger the value of k for 
which the pressure is to be evaluated (see, for instance, equation (2.17) 
and the asymptotic bcdiavior given in (2.18)), From an experimental point 
of view, thivS constitutes a considerable advantage because the amount of 
velocity measurements needed would be very limited. 
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Before leaving this section we note that in practical applications 
our analysis so far can be applied only to situations where the length of 
the cylinder is many times the diameter of the cross-section and more im- 
portant, many times the largest length scale in the spanwise direction 
(in which case the assumption of homogeneity in the x^-^iirection is valid) • 
This is usually the case at high Reynolds numb^jrs. 

2, 5 The Two-Dimensional Case , 

B^or relatively low values of the Reynolds number (Re less than 90 
in the case of the circular cylinder for instance), the large vortices 
in the wake have a correlation length in the spanwise direction many 
times the diameter of the cross-section and within a correlation length 
these vortices have straight axis parallel to the cylinder. In addition 
the random component of the velocity fluctuations is practically absent. 

In such circumstances the flow can be assumed two-dimensional and equa- 
tion (2.4) can be integrated using the two-dimensional Green's function 
G(X|Y) associated with the Laplace equation. The boundary condition at 
large distances outside the wake can be derived as in § 2.2 using (2.7) 
and the behavior of the velocity fluctuations there. The latter is given 
(upon integrating the three-dimensional behavior along the x^-direction, 
in view of the two-dimensional assumption) by: 


Note that the two-dimensional G^^e^n's function G(xjY) has a singular 
behavior for large values of 1 x-y| . In other words for very small 
values of the wavenumber k , corjjejgponding to very lar^e_^values of the 
spanwise length scale 2rr/k ^ G^x]Y;k) approaches G(x|y) only in a 
finite domain around the point Y ; this can be seen from the second re- 
lation in (2.18) and by solving (2.19) to the first order in k^ . This 
singular behavior however does not introduce any difficulty in the follow- 
ing because in (2.30) the first integral can be integrated by parts and 
the second is negligibly small as we will see in the next chapter. 
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» 

«<*c 






or 




*) *7 » 

2 2 ^ 1 *^) 
and theroforo p* decay^s like l/(xj" +x^ ) for large values of jX; 

outside and inside the woke. Hence is the solution of the following 

prob lem: 


) 

7 " 


P' 


q 




n.7p'|^ - 4n.7^v*j^ > (2.2H) 

pi jXj CO , 

Ixi 


The two-dimensional Greenes function G(XjY) is the solution of: 
7*^G - 5(X-Y) 

n.7G|(, = 0 > 


(2.29) 


7G1 -•() <3S jxl e° 


I 

G(X|Y) as defined by (2,29) is nothing but the velocity potential of the 
flow created by a negative unit source at the point Y in presence of the 
solid boundary C , The function G can therefore be constructed for a 
variety of shapes of the boundary C , using the method of images and the 
techniques of conformal mappings, as it all be *^hown in the next chapter. 

It is to be noted that the solution of (2.29) is determined up to an 
arbitrary harrr nic function, regular everywhere outside C , bounded at 
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infinity and satisfying the second condition of (2.29j. However, it can 
be easily verified that when such an arbitrary function is added to G 
the resulting expression for p' is not altered. Combining (2,28) and 
(2,29) as in § 2,3 and noting that the integral: 

J (Gn.Vp* -p'n.VG)ds(X) 

over a large circle of radius R and centered at the origin, vanishes 

i\ 

as R goes to infinity, we obtain: 


(Y,t) = JjG(xjY) q(X,t)dX + uJ G(xlY)n.V“v (X,t)d.s(X) 


(:!.30) 


The Green’s function in two dimensions behaves like log jX| for 
large values of |X; . Therefore in order to improve the convergence of 

the integral over Z in (2,30), we perform a double integration by parts, 
by using (2,5) and (2.6), in the following manner: 

-‘A 

1 J 




( 5G V 


^ij 3X.^ BX, ^ "BX.'^oX. ^ij 
11 1 1 11 


Hence : 




E 


1 J 


C+C 


R 


•J j 


Now: 


9q. . 3v. 5vl ^vl 

-1 _ J n J J 


26 


Therefore both q . . and vanish 

’ij -Xj 

on C . On G„-W , q , . is of the 

K X J 

2 

order of 1/R" and ^q. ,/dK. is of the 

i.l j 

3 

order of 1/R , consequently the inte- 
grals over vanish as R goes to 

infinity. On W , the wake width, which 


is of the order of R' 


is of the 


order of 1/R’ and 3q. ,/3X, is of the 

i.l 1 


order of 1/R 
have finally: 


3/2 



W 


and therefore the integrals over W also vanish. We 


p'(Y,t) = JJ q^.CX.t) dic + uj G(xlY)n.V^v'(X,t) ds(X) .(2.31) 


1 J 


2.6 Structure of the Source Tem Q Some Practical Considerations . 

In many practical problems an interesting and important question 
arises as to how the magnitude of the pressure fluctuations on the cy- 
linder surface, or the magnitude of the lift and drag forces varies with 
Reynolds number. A quantitative answer to such a question would require 
quantitatively detailed information about the changes which occur in the 
flow when the Reynolds number is varied. However, at this stage of our 
knowledge of the complicated flow field behind bluff bodies, a qualitative 
analysis is useful in pointing out the major factors which come into play 
in determining the pressure fluctuations on the cylinder. To this end we 
take in this section a closer look at the source term Q in the three- 
dimensional case (corresponding to high Reynolds numbers) and compare 
with the two-dimensional one. Also, we concentrate on the discrete compo- 
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nent of the flow and the corresponding surface -pressure which is dominant 
at almost all Reynolds numbers and is of considerable interest in applica- 
tion. 

First, we note that the continuity equation (the second in (2.1)) can 
be shown to be satisfied separately by the mean and the fluctuating parts 
of the velocity vector. Therefore, we have; 


BVf Bv2 


(2.32) 


and 


Sv' 

OV) 


Sv^ 


(2.33) 


Also if we denote by O)^ and u)| , i = i,2,3 the mean and fluctuating 
components of the vorticity vector o) , we have: 

^1 5x„ 5x„ 


lu: = 


Bvj^ 


9v^ 


Sx„ Bx, 


(2.34) 


a> 


Bv^ Sv^ 
3 Bx^ Bx^ 


j 


Upon taking the Fourier transform of both sides of each relation in (2.34) 
and defining (i = l,2,3) as: 


+ CO 


-ikx„ 


n 


.(X,k,t) = J co^(X,t)e ^dx^ •, i = 1,2,3 


(2.35) 


28 


we obtain; 


“ 3X., 




5V2 


:A.. 
^1.1 

^ <^X . ^X . 

1 J 


^11 ^ ^12 ^ ^ ^13 

2 Sxjix, 3 xj3=<3 


'”‘2'”'l 


" ^22 . ^^23 

- 2 Bx_Bx 

ox^ 2 3 


+ " ^32 ^ ^33 

^**2 8 x 3 ^ 

By referring to (2.6) and noting that = 0 , we have: 


2 Bx, Bx, 


0 q,, o q 

.1 ii±_ j, 

Bx Bx, 

2 1 Bx, 


^"<V3 + ''K ■ ''i''3> 


''2''3 - ''2''? 


0 *> 

+ T PCv^” - v^“) . 


Let us designate by F[a] the Fourier transform with respect to of 

any quantity Ci . We have by applying the Fourier transform to q and 
noting, in view of the assumption of homogeneity in the x^ -direction, that 
— 2 

vjv^ , v,jV^ , and v^ are independent of x^ : 


Q = 


5 Qll ^ Qj2 




-.1 ^ 




-Fik^2p(v^V3+F[v|v-]) 


n ^ 2p(v2V3+F[v'vp) 


k^oF[v^^3 , 


where 


-f- 00 


-ikx. 


Qij = J q^j(x,t)e dx3 ; i,j=i,2 


(2.37) 


or, using (2,32) ; 




3V f)V ^ 

’’ ■ S^ + 21kp(v,^+V2jj^+rt5|-(v;vp+-j|-(v^v')]) -k^F[v^-] 


The second term on the right-hand side can be written as: 
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Jv' cvi 


5v' 


2ikp(Vi +v| 


5v ^ 
^2 


or, using (2.33), (2.34), and (2.36) as; 


21kpCvj(lkVj -Li^) +V2(ikV^ + r2j) + 


^v' 




and finally as; 


2ikp { lk(VjV^ + v^V^) + -’^ 1 -^ 2 ^ -'likPCv^^] 


+ ^ ikF[v|^ + v,p] +F[v^u)| -vj^tju^]} 


Consequently Q is given by: 

- 

-' Q • • _ _ 

Q ^ + 2ikp(v^n^ +2ikoF[v^coj -v|(«p 

i” J 

- 2k"p(v^V^ +v,V^) - k^p F[vj^ d-v^^] . (2.38) 

We emphasise that the first term on the right-hand side of (2.38) is the 
Slime of four terms only and Q^. is given by: 

Q. , = p(2v,V. + f[v!v 1 - vTiT]) ; i,j = l,2 (2.39) 

In order to appreciate the expression of Q in the form (2,38), another 
feature of the discrete component of the fluctuating wake is worth ela- 
borating upon. We have said in a previous section that this component 
results from the interaction of the two shear layers emanating from the sur- 
face of the cylinder, and has a finite length scale in the spanwise direc- 
tion whose value varies with Reynolds number. While the interaction of 
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the two shear layers is a two-dimensional pheuotaenon, it is ohvimis that 
the finite spanwise length scale is the result of some three-dimensional 
effect. This effect which can be expected to be due to the presence of 
three-dimensional disturbances in the flow, takes place in each of the 
two shear layers separately. The evidence for this is provided by the 
experiment of Humphreys (1960) in the particular case of a circular cy- 
linder at Reynolds numbers close to the critical value. Humphreys ob- 
served a periodic cellular structure along the length of the cylinder 
which was not altered by the introduction of a splitter plate in the wake, 
the effect of which is normally to ])revent the interaction of the two 
shear layers (more discussion of this point is given in the next chapter). 
As tor the nature of the effect of the three-dinuMis i nna 1 disturbances on 
the individual shear layers, one might think that it is restricted to 
the generation of the small scale fluctuations which lead to the random 
component. However, since potentially amplif table disturbances with wave- 
vectors in all directions are usually present, the ones with wavevectors 
making a large angle with the direction of: the mean flov; have wavelengths 
large compared with the thickness of the shear layers. This is due to 
the fact that in the direction of such wavevectors the velocity jump 
across the layer is relatively small and correspondingly the wavelengths 
of the most amplified disturbances are likely to be large. Furthermore, 
when such disturbances grow to the extent of being able to change the 
direction o" the flow streamlines from that of the main flow, which is in 
the plane normal to the cylinder, secondary flows in the form of flow-wise 
vorticity are generated which strengthen the spanwise non-uniformity and 

Flow-wise vorticity is generated whenever a shear flow is turned in a 

plane normal to the velocity gradient. 
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create in a remarkable and puzzling manner a stable situation whose energy 
is continuously supplied by the mean shear flow. This flow-wise vorticity 
which is absent in the two-dimensional case constitutes, in addition to 
the finite character of the axial correlation length, the main distin- 
guishing features of the flow at high Reynolds numbers, To see how this 
is reflected in the source function Q we consider the integral over 
in the solution (2.15) namely: 


JJgCxIyjR) Q(X,t) dX 


Referring to (2.38) and integrating by parts the first term as in the 
two-dimensional case we obtain; 


if o^G(x'rY;k) 
1.1 3x,oX. 


Q. .(X,t)dX 


+ 2ipk jj G(xlY;k) { 

- pk^ JjG(xlY;k){2(^V^+^V2) +F[v|^+v^^]}dX . (2.40) 

The first integral in (2.40) is the same as the one in (2.31) except that 
in (2.40) only the component with wavenumber k of is taken, and 

correspondingly G(XlY;k) instead of G(XjY) is used (note G(X|Y) is 

— * 1 —4 

nothing but G(XjY;0) ). Hence varying the Reynolds number means, for 
the discrete component, varying the value of k , higher values of k 
being in general associated with higher values of Re . Therefore, this 
term which can be thought of as a quasi -two -dimensional term reflects only 
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t he el feet of the l inite spanwise correlation length of the large vor- 
tices at high Reynolds number. The second integral in (2.40) represents 
the contribution of the fluctuating flow-wise vorticity to the pressure 
fluctuations. In fact, the mean flow being steady the first term in that 
integral reflects the departure of this flow-wise vorticity Irom its mean 
position, since (v.fJj vanishes when the vorticity vector in the 

plane coincides with the mean velocity vector. The second tern 

is a nonlinear effect resulting from the non-alignment of the two-dimen- 
sional couiponents of the velocity and vorticity fluctuations vectors. 

The third integral is an additional term involving only the components of 
the velocity fluctuations in the plane (Xj,X.,) . 

in conclusion we note that the solution of the problem as given in 
(2. lb) or in the form (2.2b) and (2.27) is a function of time. In prac- 
tical problems one is interested in averaged quantities like correlations, 
intensities, and frequency spectra. Therefor^, we define in general the 
twp-point space -time correlation of the surface-pressure associated with 
the wavenumber k as: 


Rpp(Y,Y,k,r) - P(Y,k,t) P‘-XY’,k,t+r) 


(2.41) 


where the star indicates that the complex conjugate should be taken. 
Similarly for the lift and drag forces we define: 


- L(k,t) I2'(k,t4-T) 


(2.42) 


I^j)(k,r) = D(k,t) D-(k,t +r) 


When the formal expressions of P , L and D are used in (2.41) and 
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C’,42'1, we obtnln sirniK of double integrals over terms involving thefactor 
G(Xt Y;k)G(X' 1 Y' ;k) or G^jCxlYjklG^CX' 1 y' ; k) , and averaged products of 
the form: 


R^g(X,X',k,T) ^ A(X,k,t)B*fX',k,t4r) 


where A and B are the Fourier transforms of some physical quantities 

‘■4 

a(x,t) and b(x', t+r) . In accordance with the definition of Fourier 
transforms of quantities stationary with respect to the variable of 
transformation we have: 


R^j^(X,X',k,T) 


m f/2 

-• J J a (x, t)b(x ' , t +T)e 

-A/ 2 -je/2 


IkCx^-x^) 




where I is very large but finite (for instance the length of the cy 
1 inder) , or; 


in m 


ikCx^-x^ 


dXjdK;; 


1/2 

= SL I R , (X,X',§,r)e"^'^^d5 

Ur?. 


Consequently, the spectrum densities per unit length of the cylinder 
Rpp/^ , > find given as sums of double integrals over 

terms involving quantities of the form F [R^^] . Therefore, from an ex- 
perimental point of view the quantities to be measured, in order to deter- 
mine the pressure characteristics, are the various cross -correlation func- 
tions R^^ ''X,X' .’i,r) . 
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? ,7 Mmltarlons of the Incompressible Flow . 

The general analysis given in the previous sections is based on the 
incompressible form of the Navier-Stokes equations. In a steady flow it 
is well known that variations in the density of the fluid can be neglected 
if the Mach number is very small. However, for a time -dependent flow this 
condition is not sufficient, for no matter how slight the compressibility 
of the fluid is, compression waves make their appearance in the medium and 
affect a change in the incompressible flow condition. Nevertheless, it 
can be shown that within a length scale of the motion, and if that length 
scale is much smaller than the wavelength of these compression waves, the 
dynamics of the flow can be described as if the fluid were incompressible. 
More specifically if X and t are respectively a length and a time 
scale of the velocity fluctuations the condition for local incompressibility 
is : 

X« a T , C1.43) 

where is a typical speed of sound in the medium (see Landau and 

Lifshitz 1959, ?10 or Batchelor 1967, §3,6) • Since the Navier-Stokes 
equations are simply the expression of the dynamics of the flow within a 
small neighborhood of a given point the above condition (together with 
that of small Mach number) is sufficient for the validity of equations 
(2.1). However when these equations are combined together (to yield equa- 
tion (2.4) and later (2.10)) an equation of the elliptic type is obtained 
which says that the pressure fluctuations at a point are not determined 
only by the local velocity field, but by the velocity field (both mean 
and fluctuating parts) throughout the entire unsteady region of the flow, 
the effect of each point, whether close or remote, being propagated with 
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iin inl'inlte speed. For this to be valid it remains to be shown that the 
major conti ibiition to the pressure fluctuations as given by the solution 
of that elliptic equation (i.e., equations (2.15) and (2. 3D) comes from 
points within a distance of the order of the length scale X , 

In the wake behind a cylinder several length and time scales are 
present in general. However since we have considered each length scale 
^ in the axial direction separately, we can check the . Mition (2. '+31 
for each separate value of k . When k corresponds to . .e discrete 
component of the wake the time scale can be taken as the inverse of the 
shedding frequency f of the vortices and the length scale in the plane 
(Xj as the longitudinal spacing a of these vortices. Both a and 
f depend on the shape of the cylinder and for a given cylinder on the 
value of the Reynolds number. However for the sake of illustration we 
consider briefly the case of a circular cylinder. Then (2.43) can be 
written as; 


or as; 
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fmd finally as: 
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where St is the Strouhal number defined by St = ~ . The lateral spac- 

^co 

ing of the vortices is generally of the order of the diameter d and 
therefore the ratio d/a is approximately equal to 0.28 (the Karman ratio). 
On the other hand St is always between the values 0.12 and 0.21 (Roshko 
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1954"). Consequently the above condition is the same as that of small Mach 
number. As for the relative importance of the contributions to the pres- 
sure fluctuations (on the surface of the cylinder tjr instance) from 
different regions of the flow field, it can be shown that, due to the 
rapid convergence of the integrals involved, the major one comes from 
points within a distance of the order of the longitudinal spacing of the 
vortices. In fact is is shown in the next chapter that this is true even 
for the least convergent case, namely the two-dimensional one. 

When che random component is considered, for each value of k there 
exist several length and time scales in the (Xj,X 2 ) plane. Nevertheless 
a length scale of the velocity fluctuations in lov; Mach number turbulence 
is much smaller than the wavelength of the sound having the frequency of 
those fluctuations. Therefore (2.43) can be assumed to be satiffied for 
all the scales present. On the other hand for a given value of k the 
rate of convergence of the integrals in the solution is fixed since the 
Green's function G(XjY;k) is fixed* Hence for each k there exist fre- 
quencies which are high enough for their wavelengths to be comparable with 
or smaller than the size of the region over which these integrals converge 
satisfactorily. For such frequencies the incompressible model breaks down 
and the compressibility of the medium must be taken into account. 
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CHAPTER 3 


THE CIRCUrAR CYLINDER IN UNIFOM CROSS FLOW 
3>1 Introductl on > 

We shall now concern ourselves with the case of a circular cylindt‘r 
in uniform cross flow. Our present state of knowledge of the nature of 
the unsteady forces acting on the surface of a circular cylinder in uni- 
form cross flow Is duo primarily to a number of experimental investigation:; 
which were conducted and reported in the literature during the last twenty 
years. These investigations consisted in the first place of measuring 
directly the forces in question, either locally as pressure distributions 
or on a shr>rt segment of the cylinder as total lift and drag forces. In 
some instances direct measurements wore made of the moment acting on tho 
cylinder, from which values of the average forces were deduced. The ob- 
jective behind these various measurements was to provide valuable data 
necesvsary for the design of cylindrical structures subject to cross winds, 
and also to provide Information needed in the evaluation of the acoustic 
intensity radiated as aeollan tones. Concurrently, but quite independently, 
considerable effort was directed toward achieving a better understanding of 
the unsteady wake structure behind the cylinder. This effort was a contin- 
uation of an earlier interest in the subject which began after the dis- 
covery of the asymmetric arrangement of the vortices in the wake by 
Ahlborn and Benard, and the subsequent work on the stability of such ar- 
rangement by von Karman. Except in few cases, these two aspects of the 
problem, namely the wake structure on one hand and the unsteady loading 
on the cylinder on the other, continued to be treated separately and no 
attempt was ever made to study their interdependence. This is perhaps not 
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surprising because the study of the development of the vortices In the wake, 
their stability, their convection, and subsequent decay downstream was con- 
sidered of a fundamental interest, while the determination of the unsteady 
loading on the cylinder was motivated by purely practical problems. In 
fact, the part of the wake which is intimately connected to that loaduig, 
namely the near-wake region, where the vortices are first formed and which 
is sometimes called the 'formation region* was seldom investigated in de- 
tail by the students of the first aspect. The few papers which considered 
the relationship between the structure of the near-wake and the unsteady 
forces were motivated by the inconsistency of the available data on tliese 
forces: measurements made under slightly different conditions yielded re- 

sults with substantial discrepancies. Such discrepancies were suspected 
to be due, in addition to differences in the modes of measurements and the 
imperfections of the measuring devices, to differences in the free-stream 
turbulence level and the details of the experimental setup, and an attempt 
was made to explain the effect of these differences on the unsteady forces 
through their effect on the properties of the near-wake region. Although 
this attempt was only partially successful, being limited to qualitiative 
arguments, it led to substantial Insight into the mechanics of the 'for- 
mation region.' Numerical solutions of the complete Navier-Stokes equa- 
tions were equally attempted, but these cannot be expected to shed any 
light on controversial issues of the type in question since they depend on 
the manner in which the flow is perturbed and are based on purely two-di- 
mensional models. 

We present the theory of the preceding chapter as an alternative ap- 
proach to the problem which clarifies the relationship between the flow 
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structure in the uc.ir-wPkt' and the unsteady loading on the cylinder and 

puts it on a more firm and rigorous basis. This will enable us to recon- 
sider some previously published incorrect results and will lead us to the 
discussion of some interesting features of the flow in the wake. In the 
next section a description of the flow field in the near-wake is given, 
together with a review of the available data on the unsteady loading and 
the available work relating the two aspects of the problem. This is fol- 
lowed by a detailed study of a representative case of the low Reynolds 
number regime. Finally, the high Reynolds number case is discussed in 
the last section. 

3 . 2 Desc.ri lit Ion of the Flow Field - Review of Previous Work . 

The entire range of Reynolds number Re (based on the cylinder diame- 
ter), insofar as the nature of the fluctuating wake is concerned, can be 
broadly classified into three ranges (Roshko 1954). The first range ex- 
tends from a value of Re between 30 and 40 to some value between 
15i) and 200 and is characterised by a stable and regular vortex street 
extending far downstream. The fluctuating energy of the flow is concen- 
trated into discrete frequencies and decays dovmstream by the mere action 
of viscosity; no turbulence occurs in the flow. The second range which 
extends from Re between 150 and 200 to Re between 300 and 400 

is a transition range characterized by irregular bursts in the velocity 
signal from a hot-wire placed in the flow indicating a laminar-turbulence 
transition and making the determination of the frequency of the dominant 
component rather difficult. The third range called the irregular range 
extends from Re between 300 and 400 to a value near 10^ the highest 
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value of Reynolds number investigated up to date. In this range turbulent 
f 1 .ictuations accompany the periodic fluctuations; however the frequency of 

the latter is easily detected over most of the range. There is transfer 
of energy from the discrete components to the random ones as the fluctua- 
tions are convected downstream, and the wake decays by the combined action 
of viscous and turbulent stresses. It can be noted right away that the 
reason for which the limits between the various ranges are not well de- 
fined lies in the sensitivity of the flow to the experimental environ- 
ment. For instance, the value of Re which marks the first appearance 
of fluctuations in the wake depends mainly on the ratio of the tunnel 
x^;idth to the cylinder diameter, being higher for smaller values of this 
ratio. On the other hand varying the level of turbulence in the free- 
stream can vary the limits between the first and second ranges and be- 
tween the second and third ranges. We now consider in detail each of the 
three ranges separately. 

The Stable Range . 

For values of Re below 30 or 40 (and above a value near 6) the 
flow in the wake is steady and consists of two standing oppositely rota- 
ting eddies at the back of the cylinder 
and a laminar trail immediately down- 
stream as shown in the figure. The 
vorticity generated in the boundary 
layer on the forward face of the cy- 
linder goes in part into the standing 
eddies and in part into the trail form- 
ing the laminar wake (see Batchelor 1967, § 4.12). Foppl (see Goldstein 
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1938) showed, using an Idealized model of two potential flow vortices that 
the symmetrical arrangement of the two eddies is unstable for asymmetric 
disturbances. It is an experimental fact (Kovasznay 1949, Taneda 1956) 
that when a value of Re between 30 or 40 is reached the laminar 
wake shows signs of instability before the standing eddies do. Regular 
sinusoidal fluctuations antisymmetric in the longitudinal velocity develop 
and are convected downstream in the form of a vortex screet. As Re is 
increased, these fluctuations move closer to the cylinder and at a value 
of about 45 they begin to affect the tip of the standing eddies which 
begin to oscillate laterally assuming an as 3 mimetric position. Below a 
value of Re between 90 and 110 (Tritton 1959, 1971) the standing ed- 
dies, although oscillating, remain attached and fluid does not leave them; 
above that value fluid moves continually out into the street, and Tritton 
speaks of two different modes of vortex shedding, a 'low-speed mode' which 
has its origin in the instability of the laminar wake and a 'high-speed 
mode' in which shedding starts in the immediate vicinity of the cylinder. 
Tritton (1959) discovered these two modes when studying the variation of 
the Strouhal number with Reynolds number, he observed a discontinuity in 
the St vs. Re curve at Re equal to 90, separating two different laws 

of variation of St with Re , and also irregularities in the street when 

Re was close to that value. Later Caster (1969, 1971) found that such 
discontinuity or irregularities could arise from slight nonuniformities in 
the cylinder or in the free-stream in a direction along the span. Berger 

(see Wille, 1966) found another mode, t’ ' 'basic mode,' for Re above 

120 similar to the 'low-speed mode' of Tritton. Tritton (1971) argued 
that the 'basic mode' of Berger could possibly occur throughout the whole 
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stable range instead of his two modes depending on the level of distur- 
bances in the free -stream. These different findings were for some time 
the subject of a controversy which does not seem to have been settled yet 
(see the review by Berger and Wille, ^ ,0 • 

The basic mechanism involved in the development of the laminar vortex 
street is two-dimensional; however three-dimensionalities in the wake in 
the form of slantwise shedding or waviness in the vortex filaments have 
been observed in this range of Reynolds number. These observations vary 
from one experiment to another especially as far as the degree of inclina- 
tion of the vortices to the cylinder axis and the existence of waviness 
are concerned, which emphasizes once more the sensitivity of the flow to 
the experimental conditions and the level of disturbances in the free- 
stream. Gerrard (1966a) and Berger and Wille (1972) summarized the ob- 
servations of various experimenters, and the picture which emerges indi- 
cates that for values of Re not exceeding 60 the flow is stable for 
small three-dimensional disturbances or slight non-uniformities in the 
cylinder or in the flow upstream and is truly two-dimensional. This is 
supported by experiments in both water tanks and wind tunnels (Kovasznay, 
1949; Taneda,1952; and Phillips, 1956). Gerrard (1966a) further concludes, 
on the basis of his own experiment and his analysis of the occurence of 
three-dimensionality that for values of Re below that of Tritton's 
transition ’’the wake is intrinsically stable and would exhibit a stable 
two-dimensional character if the flow and model arrangement were two-di- 
mensional," This suggests that for such values of Re the attached ed- 
dies behind the cylinder provide at least in the near-wake region a sta- 
bilizing factor against three-dimensional disturbances. On the other hand, 
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for values of Re above that of the transition the shear layers which 

spring freely from the sides of the cylinder are more vulnerable to such 

disturbances and slantwise shedding or waviness in the vortex filaments 

occur inevitably as it is observed by most experimenters. 

Results on the fluctuating loads in this first range of Reynolds 

number are scarce, in particular, the ones from direct measurements. The 

only available data are due to Tanida et al. (1973), and their Figure 

3(b) is reproduced here for reference. In this figure |c' | is the 

L 



amplitude of the fluctuating lift coefficient defined by: 

where L is the lift per unit length. Tanida et al. made their mea- 
surements in an oil tank (for the low values of Re ) on a 100 mm long 
central segment of a cylinder 352 mm long and 30 mm in diameter. The 
main feature of their results is the low level of lift fluctuations at low 
Reynolds numbers. Phillips (1956) using the data of Kovasznay (1949) at 
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Re equal Lo 56 , calculated the amplitude of the flucLuatiuR lift and 
drag and obtained the rather high value of iCj| = 0.76 for the lift. 
Jordan and Fromm (1972) solved numerically the complete time -dependent 
incompressible Navier-Stokes equations in two dimensions and obtained at 
Re equal to 100 the value Ic, ] =0.27 . These variations in the results 
will be discussed at length later in this chapter. 

The Transition Range . 

When Re reaches a value between 150 and 200 , the flow in the 

near-wake, though still laminar, contains low-frequency irregularities 
which become more violent downstream and eventually render the far-wake 
turbulent (Bloor 1964). These low-frequency irregularities are believed 
to be due to three-dimensionalities in the flow. However, small bursts 
of turbulence do not make their appearance be lore a value ol Re of 300 
is reached (Bloor 1964) . Thesc^ occur at random in time and in a direction 
parallel to the cylinder axis making the vortex filaments in the street 
partly laminar, partly turbulent until a value of Re equal to 400 is 
reached. It has been argued by Roshko (1954) and Roshko and Fiszdon 
(1967) and observed by Bloor (1964) that these turbulent bursts occur 
right before the end of the formation region which is defined by Bloor 
as the beginning of the (turbulent) vortex street and by Roshko and Fisz- 
don as being located near the closure point or the tip of the two standing 
eddies of the corresponding mean flow. As a result of the low-frequency 
irregularities associated with large scale three-dimensionalities and the 
intermittent turbulent bursts, the flow in the near-wake in this transi- 

"imai (1964) has shown that up to a value of Reynolds number of at least 
6000 the mean flow pattern is similar to that at low Re where two 
standing eddies exist right behind the cylinder. 
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tioti range is obviously three-dimensional. No information is available on 
the nature of the unsteady loading on the cylinder In this range except 
the numerical results of Jordan and Fromm (1972) at Re = 400 . They ob- 
tained the value 0,75 for the amplitude of the lift coefficient. 

The Irregular Range . 

When a value of Re of about 400 is exceeded the turbulent bursts 
start to occur more systematically before the end of the formation region 
so all vortices in the street are turbulent on formation (Bloor 1964) . 

The low-frequency irregularities continue to exist in the formation region 
and farther downstream with substantial reduction in intensity in the region 

in between. This state of affairs continues until a value of Re of 
3 

about 1.3X10 is reached when sinusoidal waves of the Tollmien-Schllch- 

ting type begin to precede the turbulent bursts (Bloor 1964). These, 

called by Bloor, "transition waves " have a definite frequency f^ which 

is greater than the fundamental shedding frequency and varies (for a con- 

3/2 

stant cylinder diameter) as ' . The ratio of f. to the fundamental 

00 i- 

1 ,. 

frequency is proportional to Re^ at high Reynolds number and it is equal 

3 3 3 

to 2.5 at Re = 1.3X10 and to 8 at Re = 5X10 . Below 1.3X10 the 

transition waves are not detected and as suggested by Roshko and Fiszdon 

(1967) it is possible that, since their frequency is very close to the 

shedding frequency in this range, there is a coupling between the two 

3 

modes, in particular when Re is close to ^lOO . Above 1.3X10 the 

low-frequency irregularities continue to exist in the formation region but 

3 

disappear in the wake downstream. As Pe is increased from 1.3X10 to 
3 

8X10 (Bloor 1964) the transition points in the separated shear layers 

move upstream rendering the downstream part of these shear layers more and 

3 

more turbulent. Beyond Re = 8X10 the transition waves are no longer visible 
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inside the v?ake, the laminar portion of the separated shear layers being 

rapidly followed by turbulence (Bloor 1964). Transition to turbulence 

continues to move upstream along the shear layers as Re is hicreased 

until it becomes very close to the separation points on the surface of 

5 

the cylinder. When a value of Re of about 2xlo , the so-called cri- 
tical value, is reached and while the separation is still laminar, 
transition is followed by reattachraent which in turn is followed 
by turbulent separation on the back of the cylinder. The wake is 
narrower, there is a sudden fall in the drag coefficient and a loss of 
periodicity in the wake. This continues until the transcritical regime is 
reached (Re>3,5X10 ) in which transition to turbulence precedes separa- 
tion (Roshko 1961) and where a definte vortex shedding occurs as at sub- 

critical Reynolds numbers. The critical or supercritical regime 
5 6 

(2X10 < Re <3.5X10 ) is characterized by a great sensitivity to surface 

roughness and free-stream turbulence, one manifestation of which is the 

loss of regular vortex shedding caused by gross non -uniformity along the 

length of the cylinder. Using highly polished cylinders Bearman (1969) 

5 5 

conducted interesting experiments in the range 10 <Re< 7.5X10 and 
found that by carefully and frequently cleaning the cylinder during the 
experiment regular vortex shedding could continue up to a Reynolds number 
of 5.5X10^ . 

There is less disagreement among the various experimenters about the 
existence and nature of three-dimensionalities in the irregular range than 
in the stable or transition range. As noted by Gerrard (1966a), '^this is 
perhaps because the chaotic nature of the three-dimensional structure al- 
lows less precise description of the phenomena." Most workers describe 
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the three-dimensionality In terms of a spanwise correlation length, i.e., 
the distance in the axial direction over which the fluctuations of some 
property of the flow are well correlated. By correlating velocity fluc- 
tuations Roshko (1954) finds a value of 3 diameters at Re = 500 , 
Prendergast (1958) correlated surface-pressure fluctuations but his corre- 
lation functions did not tend to zero for large separations. El Baroudl 
(1960) repeated the same experiment but instead correlated velocity fluc- 
tuations near the shoulder of the cylinder and found a correlation length 

4 4 

which increases slowly with Reynolds numbe* in the range 10 -'Re <4.5/10 , 

bf 

being about 3 diameters at Re = 10 and about 6 diameters at Re-4,5xlO , 

3 

Phillips (1956) gives a value of 3 diameters at Re ~ 5x10 . Mattingly 

(1962) observed a spanwise periodicity on the surface of the cylinder in 
4 5 

the range 10 <Re<10 with a v^avelength of a few diameters. Using 
threads attached to the cylind:;r Humphreys (1960) observed a periodic 
cellular pattern with a wavelength of about 1.5 diameters at the critical 
Reynolds number. Gerrard (1966a) describes the nature of three-dimension- 
ality in the irregular range as a combination of randomness associated 
with the small-scale turbulent structure and a ’’more gentle variation of 
the low Reynolds number type” probably associated with the low-frequency 
irregularities . 

The bulk of the experimental data available concerning the unsteady 

loading on the cylinder corresponds to values of Re in the irregular 

3 

range, and more specifically to values of Re greater than 2X10 . Yet 

from these data it is not possible to draw any firm conclusion about the 
”true” level of the loading or its real behavior when the Reynolds number 
is varied across the range. This is perhaps best illustrated in Figure 15 


49 


J.0| 


ai^ 


0 . 6 ] 


0.4i 


0.2 


» r » T » f T T T 
|G:(?8i CHollMing 
max tmpi 


CiUW 

/Vrm!"* ' 


K: niol rms 
Discs, 3.1 


/ 


Hr [31! 
Vents sealad 


Kt 11101 rms 
Discs, IS 1 
or vents closed 


f rms miKnenj^ 



\ C;I28Jr>5KJ 

■ max amol 


/ , ' H: [311 

/ vents, rms moment 

/ 

^ y' K: [1101 rms 

Ho discs; vents open 


I iir 10 ^ 10 ^ 

Re 


' \ /- m« ampl| 

)'- 

^ F: IMI rmj 

1.74:1 SKliotv, gips I 


FIG, 15. OSCILLATORY LIFT COEFFICIENT. 


of Morkovin's (1964) review which is reproduced here and in which the root 
mean square value or the maximum amplitude of the unsteady lift coefficient 
is shown. In this figure the results of Gerrard (1961), Keefe (1961), 
Humphreys (1960), Goldman (1958), and Fung (1960) are summarized and in- 
dicated symbolically by G[lll] , K[110] , H[3l], G[28], and F[30] respectively. 
Gerrard measured the r.m.s, value of the local presare fluctuations at the 
shedding frequency and assuming a constant phase on each half (upper and 
lower) of the cylinder surface integrated to obtain the lift. Keefe mea- 
sured the forces on a central segment of 1 diameter in length and examined 
the effects of fastening to the cylinder ends two circular discs (in order 
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to redliro tho intc*rti*rc‘no(* of the w<ill boundary layers) of dlametc-r five 
times that of the cylinder and separated by a distance of ’i or IS cy- 
linder diameters and also the effect of sealing the clearance holes at 
biUh cylind(‘r ends. Hiimphreys* results are based on measurements of the 
moment on a cant ilevc^red cylinder; he also studied the effect of changing 
the end conditions. Fung measured the fluctuating forces on a short seg- 
ment of 1 , 7'f diamet t*rs at supercri tical Reynolds numbers and noted the 
effect of the gaps betw^uui the foree transducer and the rest of the cy- 
linder, but did not recognize the effect of the length of his transducer 
being larger than the correlation length at these high Reynolds numbers. 

To these various data sets, others can he added but only to complicate 
further tlu» general picture, the rtsKson being, as we have stated earlier, 
the stmsitivity of pressure or forces measuxisaents to the dc^ tails of tlu^ 
experimental setup, the Imperfections of the measuring devices and the 
level of disturbances in the free -stream. Results not mentioned above in- 
clude tliose of Macovsky (19^38), McGregor <19^37; who measured local pressure 
fluctuations, Jones C‘t al (1969) at supercritical and transcritica 1 Rey- 
nolds numbers, and also at these high Reynolds numbers those of Schmidt 
(1965, 1966). These latter (i.e., the high Re data, in particular those 
at supercritical Re ) shoxc in addition a strong dependence on the state 
of cleanliness of the cylinder which is a consequence of the extreme sen- 
sitivity of the flow to local surface roughness. 

Among the various factors which can have a major impact on the mea- 
sured values of the unsteady forces the free-stream turbulence is the only 
one which received some attention by the workers in the field. Gerrard 
(1965) puzzled by the large discrepancies between his results (of 1961) 
and those of Keefe (1961) and others, and by the low level of lift fluctua- 
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flons he observed at Re “ 4x10 , investigated the effect of increasing 

the turbulence in his wind tunnel. Instead of measuring again the pressure 
fluctuations he measured the velocity fluctuations near the shoulder of 
the cylinder and found an increase (when the wind tunnel turbulence was 
increased) in the intensity of these velocity fluctuations similar to the 
discrepancy between his results and those of Keefe, indicating that this 
discrepancy is perhaps due to differences in the level of disturbances in 
the free-stream. He also recognized the role of characteristic length 
played by the length of the fomaLion region but stopped short of consi- 
dering the effect of an increase in free-stream turbulence to be equiva- 
lent to an effective increase in because the shedding frequency re- 

mained unchanged. Later he (Gerrard 1966b) discovered a second character- 
istic length, the length to which the shear layers diffuse at the end of 
the formation region and as a consequence achieved a substantial insight 
into the mechanics of the formation of vortices in the wake. On the basis 
of the two characteristic lengths Gerrard could explain why the Strouhal 
number remained constant over a wide range of Reynolds numbers and why 
the shedding frequency was insensitive to changes in the free-stream tur- 
bulence level. However, on the basis of the same arguments Gerrard 

could not explain the great change he observed (but Keefe did not) in the 

3 4 

level of fluctuating lift when Re was varied from 4X10 to 7X10 
This led him to suggest that possibly at the lower Reynolds number and in 
the absence of free-stream turbulence the two separated shear layers de- 
velop independently of each other creating a symmetrical formation region, 
the periodic wake downstream being formed in a manner similar to that at 
Reynolds numbers below 90 . The idea of a symmetrical formation region 
was, in fact, inspired by the results of his potential flow model (Gerrard 


52 


1967a) in which the shear layers separating from the sides of the eylindtr 
were modelled by a distribution of elementary vortices moving under the 
action of the flow past the body and the velocity field of the vortices. 
This model however, when forced to generate a fluctuating lift in agree- 
ment with Gerrard's data at low Re , gave a formation region much larger 
than the value measured by Bloor (1964). Some comments on the model in 
question will be given at the end of the present chapter as part of our 
general discussion. 

Finally, mention should be ma'’ of the few data points obtained by 

3 4 

Tanida et al. (1973) in the range 2V10 <Re<10 (using the same test 
cylinder as at low Re but with water as the medium) which are shown in 
the figure already introduced in connection with the low Reynolds number 
range, and also the value 0.95 for the amplitude of the lift coefficient 
found in the computations of Jordan and Fromm (1972) at Re = 10^ . 

3.3 A Representative Case of the Low-Speed Mode of Vortex Shedding . 

In this section we consider in detail the problem of determining ana- 
lytically the unsteady pressure distribution on the surface of the cy- 
linder when the Reynolds number Re is equal to 56 . At this Reynolds 
number Kovasznay (1949) investigated in some detail the velocity field in 
the wake and his data will be used to generate the essential terms needed 
for the application of the theory of the preceding chapter. Also at this 
Re the wake structure and the nature of the relationship between the near- 
wake and the surface -pressure are typical of the whole range below Re = 90 , 
namely the vortex street develops as a result of the instability of the 
laminar wake and is stable and laminar. The velocity fluctuations are 
negligibly small near the cylinder and reach a maximum intensity some 
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distance downstream (Kovasznay 1949). The flow is two-dimensional and 
therefore the results of 2.5 are applicable. In addition, in this range 
of Re we. can make the assumption: 

I 2“^ 

n,7p'l =un.7v' 

|c 

This is based on the following argument: since the vortex street develops 

some distance downstream of the cylinder, the unsteady part cu' of the 
vorticity vector vanishes in the neighborhood of the cylinder, the small 
velocity fluctuations there being irrota tional and induced by the vortical 
street further downstream; nevertheless at the boundary C and as a re- 
sult of the no-slip condition additional unsteady vorticity might be gen- 
erated locally, in other words, while the velocity fluctuations themselves 
are vanishingly small there their gradient can be large; this latter pos- 
sibility however can be eliminated on the assumption that the frequency of 
the oscillatory motion in the wake is high, i.e., the "relaxation time" 
associated with the viscous flow is of the same order as, or even greater 
than, the time scale of the oscillations, and therefore the viscous ef- 
fects do not have the time to respond to the full no-slip condition. 

From this it follows that co* vanishes at and near the boundary C , 
and since 



2 ^. 
7 v' 


= -Curlo)' 


This is similar to the argument sometimes invoked in connection with the 
Kutta condition at the trailing edge of a thin airfoil shedding vortices 
at high frequencies. Previous experience has shown that relaxation of 
the Kutta condition in such circumstances yields results in better agree- 
ment with observations (see for instance Davis 1974). The question re- 
mains however whether the above condition is satisfied for values of Re 
above 90. It is interesting to note that in all the previous_^studies 
made on the case of a turbulent boundary layer the condition n.Vp^=0 
at the solid boundary has been adopted. 
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the above condition can be considered satisfied. The resulting expression 
for the pressure p' is: 


p'(Y,t) 


rr 9 ^g(x|y) 

.IJ dX.SX, 
r ^ J 




(X,t) dX 


( 3 . 1 ) 


3 . 3.1 The Green's Function G(X,Y) . 

As previously noted, the two-dimensional Green's function G(X,Y) is 
nothing but the velocity potential of a negative unit source at Y in the 
presence of the boundary C which in the 
present case is a circle of diameter d . 

G is therefore the sum of the potential 
of that source in absence of C , plus 
that of a similar source at the inverse 
point Y' of Y with respect to C , 
and plus the potential of a source of op- 
posite sign at the origin. In other words, G can be. written as: 

G(x1y) = - ^ log |X-y'| - i log IX-Y' i + ^ log |X| , 



or as: 


2tt 


% 


2tt 


G(x1y) . - ■^logt(X,-Yp^ + (X2-Y2)^}'- ilog{(Xj-Yj/+(Xj-Y^)^)‘ 


where 




(ii/2)^Yj 

Y ^+Y ^ 

1 2 


^2 


(d/2)^Y2 

Y ^+Y ^ 

1 2 
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Since in (3.1) the derivatives of G appear rather than G itself, we 


differentiate with respect to Xj and to obtain the following expres- 

sions : 


5X^“ [(X^-Yp2+(X2-Yp2]^ (X^^+X,“) 


1 . ^'' 2 -'' 2 ^ 


(X2-Y^)-(X^-Yp 


,^2 


2 2 

S.-S 


7-7} 


o 

d“G 


2rr 


1 ^ ^(Xi-Yj)(X2-Y2) 


+ 


2(X^-Y|)(X2-Yp 


2 X 1 X 2 


[ (Xj^-Y^)-+(X 2-Y2)2]^ C(X^-Yp2+(x^_Y^)2]'^ (X 22 +X 22 ) 




and since G is a harmonic function of and we have; 


0 

c‘'G 


9 

d^G 






2 for X^ 7^ Y^ , X 2 Y 2 


When we evaluate the surface -pressure fluctuations as given by (3.1), Y 
goes onto the circle C and Yj , Y^ become respectively Y^ and Y 2 ; 
the above expressions reduce to: 




5X 


1 r (Xj-Y.) -(X.-Y )2 X^-X^ 

7=-M — — j - 7 — ^ ’ Y one (3.2) 

^ ^ [(X^-Y^)2+(X2-Y2)2]^ (X^2+x^2)2 


5^G 1 . 2(X2-Y^)(X2-Y2) 

BX,5X^ rr ^ 


2 X^X 2 


^1^^2 ^ C(X^-Y^)2+(X2-Y2)2]^ ^ (x^2+x^2)2 


Y on C 


(3.3) 


3.3,2 The Velocity Field in the Wake . 

Both the mean and the fluctuating parts of the velocity field in the 
wake at Re = 56 were investigated by Kovasznay using the hot-wire tech- 
nique. Measurements of the total mean velocity distribution and that of 


the r.m.s. value of the velocity fluctuations were made at several down- 

ORIGIN'AI- p- 


stream stations. In addition, it was observed that the fluctuating vc* In- 
city had two components with two different frequencies, one is the shedding 
frequency f , the other is the double frequency 2f , and the phases of these 
two components as functions of were determined. At the middle of 

the wake only the double frequency component was present, but this de- 
creased rapidly with distance from the wake centerline and disappeared al- 
most completely before the point of maximum fluctuations was reached. 

In order to determine the various q.. terms in (3.1) the two compo- 

ij 

nents and v^ of the mean velocity vector are needed. However, by 

examining the mean streamline pattern plotted in Kovasznay^s work and a 
similar plot made during the present work it was judged that the v^ 
component is everywhere a small fraction (perhaps lOZ or less) of the total 
mean value except in a central narrow region of width 1.5d bounded down- 
stream by Xj = 5d ; in that region the total mean velocity is a small 
fraction of the free-stream velocity ; consequently it will be assumed 
here that t erywhere v^ is an order of magnitude smaller than . As 


for V- it is to a good approximation equal to the measured value 
f 1 
— _ 2 ~ 2 

V = "^2 ^ except in the narrow region mentioned above. This region 

contains the two standing eddies in which there is a reverse flow (see 


Taneda's 1956 picture at Re = 57.7 ). Kovasznay measured v at 


Xj^ == 2d , 3 . 5d , 5d , 8d , . . . At X^ = 2d and 3 .5d the measured values on 


the wake centerline which are exactly those of v^^ , were taken here as 
negative and an interpolation (numerical, using spline functions) was made 
between these values and values off the centerline which were judged to be 
a good approximation of v^ . In this manner the distribution of v^ in 
the region of the standing eddies where v^ can be larger than v^ and 
in particular where v^ changes sign was determined. Again it is believed 
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that Vj as determinod by the data and the above procedure is accurate up 
to 107 . 

In addition to the mean velocity components the determination oi’ 

requires a knowledge of the fluctuating velocity components v^* and 
V,,* . Since the hot-wire is more sensitive to the fluctuating velocity 
component in the direction of the mean flow and since the direction of tlu^ 
latter is, to a good approximation, in the X^-direction everywhere e:crpt 
in the neighborhood of the standing eddies, the r.m.s, values as measured 
by Kovasznay can be taken as those of v^‘ except in that region; there 
a similar procedure to the one used above was adopted. As suggested by 
Kovasznay, we can write: 


,X./)cos 2TTL0px^)-ft]+cp2(X, ,X^)cos4TT[S,^(Xp-fL] 

(3.4) 


where and are respectively the amplitudes of the fundamental 

component and its first harmonic and 0^ and 9^ their phases. In 

addition, is an odd function of X^ and cp^ is an even function of 

X^ which can be assumed to vanish parabolically in the transversal direc- 
tion, i .e , , 


cp^CX^.X^) = CP2(X^,0)[1- 



(3.5) 


where X^^ is the value of X^ for which the total fluctuations are a 
maximum. Using (3.4) and (3.5) and noting that: 


In fact, Kovasznay had a -t sign inside the brackets in the expression 
corresponding to (3,4), but since the vortices are convected in the posi- 
tive X^-direction and 9^ and 0^ are increasing functions of X^ a 
- sign is more appropriate. 
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■> 1 •) ■) 
“ = + 


> 


and can be calculated from the measured r.m.s. values. These 

latter are vanishingly small for < 2d , peak at about = 7d and 
slightly off the wake centerline, and are everywhere an order of magnitude 
smaller than the mean velocity. As for V 2 ' , which was not measured by 
Kovasznay it will be determined here using the continuity equation. First 
we note that like Vj ' , v,^ ' has two components with two different fre- 
quencies. One is f . the other is 2f . It can be easily seen that the 
continuity equation must be satisfied separately by each frequency compo- 
nent; in other words, if wo denote by v^'^ and v ^'2 the components of 
Vj ' with frequencies f and 2f respectively, and by V 2 '^ and 

those of V ' , we have: 

9 


3v 


n, 

1 


+ 





0 


and 


!ll2 . 


0 


and vanish at X^ = + , in addition 

This can be seen by examining the velocity field of 
street where it is found that v,> on the wake axis 
frequency. This is also observed experimentally in 
bodies shedding vortices (see for instance Campbell 
integrating the above equations yields: 


vanishes at X^ = 0 , 
an idealized vortex 
is of pure fundamental 
the wake of cylindrical 
1957) . Consequently, 


'21 


^>x: 


v|i(Xi,e2,t)d?2 


59 


therefore : 




= -j 


SX ■ J At '^ 12 ^^!’ •■ 2 ’*'^‘^•'0 ’ 


where 


•1 ■*■/ vi2^^i’52.t)d5 


(3.6) 


wCXpX^.t) = Jii(Xi,X2)cos2rrC0^(X^).ft]+02(X^,X2)cos4n[02(Xp.ft] 


where 


(3.7) 


X^ ) 

=/ ^ 2 ^^ 1 ’^ 2^‘^"2 I 

0 

In order to evaluate inside th. „ake, values of <p^ outside the „ake 

are needed. The available data from which cp^ »as calculated cover the 
wake region up to 5 or 6 diameters from the centerline, at these bor- 
derlines the r.m.s. values of the fluctuations are very small and so is 
1>1 . However, in order to account for the remaining part of the integral, 
l.e., the part beyond . 5 or 6d . was assumed to decay like 
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l/X^” in the +X 2 -direction (see 5 2.5). 


3.3.3 The Fluctuatlns .Surface -Pressure Field - Lift and Drag Fluctuations . 

The fluctuating surface -pressure field is given by (3,1) where Y 
is taken on the circle C . By expanding (3.1) we obtain: 


= fj ^ qjjdXjdX^H- JJ 


j, SXj 


BXjSXj- 




‘2 1 




dX^ 


We also have: 


q..(X,t) = p(2v. v' + vl v! - viv' ) 

1 .1 1 .1 X j 


Therefore: 


p'(^%t) = 2p 


r ;7T''i''i“i'“^2+2o, 
X ^^1 1 


r 5 "g — , 

J ax ^ ax ^ ^1 ^2 *^^1 ^^^2 


2 

JJ aMr^n ^>'1 '*^2 - 2pJI “i “2 

^12 


P If ^ \ II BX.9X ^''1''2 ■ ''r2'""l""2 

X ^^1 S ^ 2 


(v'v' -v'v')dX,dX. 


’’ 2 . 2. 




X ^^1 


(3.9) 


In evaluating p' only the first two terms in (3.9) will be retained, 
these are thought to contain the major contributions to the real value of 
p' - The contribution from the third and fourth terms which contain the 
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mean velocity component V 2 can be expected not to exceed 10 of the 
contribution from the first two, and is perhaps of the same order of mag- 
nitude as the error involved in evaluating the first two integrals with 
Vj as determined in the last section. The contribution from the remain- 
ing terms, which are quadratic in the velocity fluctuations, is discussed 
in the Appendix. With the understanding that we are computing only an 
approximate value we write: 


p'(Y,t) = 2p 


2 

rr^j^7;-v^dx,dx, .(3.111) 


or 


p'(Y,t) = 2p 


'’"G — , „ rr — ow 


^ 'X, 

1 


dX^ dX^ 


T. 


•} O 

I ^^1 E 1 ^ “ 


+ 2o 


7 

t' i. 


*^.2 

Vj w dx^ dx^ + 2o ’ ' 








The second integral in the above expression can be integrated partially 


with respect to between some value X^^^^ of X^ such that 

^<XjQ<2d and Xj^ = + ® . At X^^ '• 0 ^ and are vanishingly small 

and since v, and are finite there, the product v, w 


1 




3Xj3X, -1 IXjj, 


is negligible. Furthermore, as X-“*+co y and w remain finite, how- 

-s2 I i 

n G 1 

decays like — - (see equation (3.3)) and the product 
^1 


ever . 




vanishes again. Therefore the integral in question vanishes and we have: 
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Vj w dXj dX., 


-* f’f’ S'^G 

p'(Y,t) = 2 d ! ' -4 V, v!dX dX +2ojj —4- 

E '^1 I '^1 '^2 


+ =’’JJ ^0x7 ST “'=’‘1 “2 

^ 1 fc. i 


(3.11 I 


The above step, i,e. the partial integration with respect to Xj is vital 
for maintaining a good accuracy of the numerical computations. The data 
which will be used to evaluate the physical functions in the above inte- 
grals cover only the part of the wake between X^, = 2d and X^^ » 40d . 
Although and i!)^ are small for X^ < 2d , their derivatives with 

respect to X^ there cannot be accounted for numerically with a good 
accuracy. In addition, the form assumed for and the numerical repre- 

sentation of '■pj are approximate; therefore while their integrated values 
and are still a good approximation, the derivatives of the latter, 

i.e., (taken numerically) are not. Finally, 

while in (3,11) only one numerical differentiation is required, in (3, 10) 
four are necessary to evaluate . 

We now define the following functions: 


0^(Xi,X2) = 


o 9 
V V 

''2 "^1 

2 2 ^ 
(X/+X2 ) 


/iCx^.x^) = 


2X2(3Xj^^-X2^) 

2 2 ^ 
(X^^+X^^) 


(3.12) 


y(x^,X2) 


2X1X2 


2 2 

(X^ +X 2 ) 


Numerical integration is always a smoothing process, 
differentiation which reduces the accuracy. 

63 


unlike numerical 


It follows from (3.2) and (3,3) that when Y lies on C : 


- - i Ca(Xj-Yj,X2-Y.,)-ia(Xj,X2)] 


2 


(3.13) 


(3.14) 


and also after one additional differentiation with respect to Xj : 


^ - k'/3(X -Y ,X,-Y^)-|g(X..X )} . (3.15) 

5>X^ 5X, i 1 i ^ 

We note that a is an even function of X^ and 8 and v are odd func- 
tions of X., . In addition, and Svj^/3Xj^ are even functions of X., , 

and are odd functions of X^ , and finally cp^ and are even 

functions of X^ . Using the equations (3.4), (3.7), (3.13), (3.14), and 
(3.15), relation (3.11) can be written in the form: 


p'(Y,t) = - 


Is. 

TT 


JJ{a(Xj-Y^,X 2 -Y 2 )--|oi(X^,X 2 )}v^Ccp;L cos 2TT(0^-ft) + 


E 


+ fp, cos 4TT(02-ft) } dXj^ dX^ 

- ^ JJ{8(Xj-Yj,X2-Y2)--|0(X^,X2)}^[0j^cos 2n(0^-ft) + 

E 

+ \l >2 cos 4TT(0^-ft)} dX^ dX^ 
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+ »X 2 -Y 2 )-|-y(Xj.X 2 )} ^ ftl)j cos 2TT(e 


+ cos4'T(e2-ft)}dXj dX^ 


1( Is oasily seen fliat after expanding the above expression, the 
which contain the product ^ Vj and vanish 
remains can be written in the form: 


p ' (Y, t) = - JJ Otj Vj Tj cos 2rr(0j -ft) dXj dX^ 
X2^ 0 


" ^ J/ V, 'P,, eos /fTTfQ., -ft) dX dX,, 


X2-0 


fr 


+ ■'if' * I 4TT(0„-ft) dX, dX„ 

TT «• 12 2 12 

S. ' “ 


■ JJ ^1 2TT(0^-ft) dXj dX, 

X2-0 


^ II ^2 ^2 ^■^(®9~ft) dXj dX2 

X2- 0 


X2^0 


i£ 

TT t) 


X2^^ 0 


Vi ^ iij cos 2rr(0j -ft) dX^ dX 2 


-ft) + 


Integra 1 s 
and what 
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where: 


^ II '"2 ^ '''2 *^^2 


X2 = 0 


-^JJ y5:«i^2COs4Tr(02-ft)dX^dX, 


>- 1 


VO 


0!^(X;Y) = 0!(X^-Y^,X2-Y2)-0!(X^-Y^,X2+Y2) 


a^CXjY) = a(X^-Y^,X2-Y2)+0!(Xj-Y^,X2+Y2) 


j 3^(X;Y) = ^(X^-Y^.X^-Y^)- ^(X^-Y^.X^+Y^) 


^^(X-.Y) = ^(X^-Y^.X^-Y^j+flCX^-Y^.X^+Y^) 


2 


(3.16) 


(3.17) 


V^(X;Y) = v(X^-Y^,X2-Y2)-V(X^-Y^,X2+Y2) 

y^aiY) = v(X^-Y^,X2-Y2) + y(X^-Y^,X2+Y2> 




(3.18) 


Finally, if v?e denote by (ju the angular frequency 2TTf and we introduce 
the non-dimensional quantitites •^/U , co, /U , cp„/U , )1), /U d , lI)„/U d , 

a(v^/uj X X 1 o 1 » 2 « 1 » 2 « 

a(x /d) • T ’ T “ 


p'(Y,t) 


2 

— cos 

rr 


X, X„ 


(*P 9 ^1 'Pi ^1 *^1 „ “1 “2 

00 CO 00 ® 1 


X2^0 


2 . rr 2 ^1 '<’1 ^1 *^'i . „ ^1 ^2 

- ^sxn OJt JJ + U^U^d " 3 (X^/d) d d ^ 




00 00 00 CO 
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dVi(X;Y) = V^(X/d;Y/d) 

and that similar identities hold for ^2 ’ ^2 ’ '^1 ^ , y . 

As is clear from (3,10) the surface-pressure fluctuations contain two 
components with two different frequencies, one is the fundamental or shed- 
ding frequency and is represented by the first two terms in (3,19); the 
other is the double frequency and is represented by the remaining four 
terms. In addition, we note that , and are odd functions 

of while ^2 ’ ^2 ’ "^2 even functions of . This implies 

that the fundamental frequency component is antisymmetric with respect to 
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to the direction of the free-stream, while the double frequency component 
is symmetric with respect to that direction. In terms of the lift and 
drag fluctuations this means that the lift fluctuates with the shedding 
frequency and the drag fluctuates with twice that frequency! 

The data of Kovasznay were used to determine the integrands ?‘n (3.19) 
as was indicated in the last section. These data were first smoothed and 
interpolations were made using spline functions; the results of the inter- 
polations were checked and proved to be excellent. Then the various func- 
— ^2 ^2 

tions (5^ , 0, , v^/U^ , - « iT . ^ ’ u“ determined in 

the domain 2 ^ X^/d ^ 20 , O^X^/d^-S using Simpson's rule of integration 

for the last two, and the aforementioned interpolating functions (which 

are twice continuously differentiable) for <> '• ' . Finally the coef- 

a{Xj/d) 

ficient of cos 0)t , sinu't, cos 2cjot and sin 2o)t (i*e, the integrals) in 
(3.19) were evaluated at 19 different points (separated by a uniform angu- 
lar distance of 10^) on the upper side of the circle C . The integration 

^1 ^2 - 

was dotxe over the domain 2^— ^20 , 0^ — --5 , using Simpson's rule and 
a mesh size of 0,5 d , This domain covers in the -direction two wave- 
lengths of the vortex street, and the major part of the unsteady region 
in the X*-direction, In the next section it will be shown that the major 
contribution to the complete integrals comes from this limited domain. 

For convenience and in order to determine the phase variations of 
the surf ace -pressure fluctuations we write (3.19) in the following form: 


pU 


2 cos 2rr(6^-f t)^?^ cos 4Tr(6^-ft) 


(3.20) 


where ’ ^2 ’ ^1 ’ ^2 functions of the angular position on the cir- 

cle C , In addition, P^ has an antisymmetric distribution with respect 
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to the X--axis, and , 6 , and 6 have symmetric distributions with 

respect to that axis. These functions were determined and are plotted In 
Figures 1 and 2. In Figure 2 the phases are plotted in such a way (by a 
translation of the phase ordinate, different in each case) that they are 
zero at the back stagnation point. It follows from the form of the phase 
distributions that as the vortices are shed in the downstream direction 
the pressure pattern at each frequency travels along the surface of the 
cylinder in the upstream direction. Note also that the intensity of the 
shedding frequency component is maximum at about 120^ from the front stag- 
nation point while that of the double frequency component is maximum on 
the back of the cylinder. These results will be discussed in a later 
section. 

The results of the integration in (3.19) can be used to determine the 
magnitude of the lift and drag fluctuations on the cylinders. If L(t) 
and D(t) designate respectively the 
lift and the fluctuating part of the 
drag per unit span we have: 

2tt 


J 

0 


p ' sin 


2tt 


9d0 


D(t) = “ 2 "/ p'cos0d0 
0 



where in the first integral only the antisymmetric part of p' has a net 
contribution while in the second only the symmetric cup does. 
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0.04 
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Effecting the integration yields: 

pU^^d 

L(t) 0.0186 -y- cos (Ut , (3.21) 

pU^^^d 

D(t) 0.0007 — 2 — ’ (3.22) 

the time origin being not the same in both expressions. From (3.21) and 
(3.22) it is seen that the magnitude of the lift fluctuations is about 
25 times larger than that of the drag fluctuations. 


3.3.4 Estimation of the Error Involved in Integrating Over a Finite Domain. 

In this section the error involved in evaluating the right-hand side 

^1 ^2 

of (3.19) by integration over the limited domain 2 ^-^^20 , is 

examined. The estimation procedure is instructive and elucidates the man- 
ner in ^^ 7 hich the velocity fluctuations in the wake induce the fluctuating 
surface -pressure field . 

Referring back to the expression of p‘ as given by (3.10) and to 
the relation (3.2) and (3.3) we note that for a fixed value of , we 
have: 




CO 


) 


a^G 1 


for 



CO 


and similarly for a fixed value of ; 



for IxJ 


CO 
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re 1 




for |Xj 1 -•* “ , 


^2 

but for a fixed ratio nf ~ we have: 


JLL ^ j:- for lx, 

5X,2 x; 


1 1 y 

^ ^ x^2 1 1 


This is due to the fact that 


is of the order of 


1 


fJX^eiX^ 


X„ 


admits a local, extremum whose, value 


and whose X -coordinate increases linearly with 


However for a fixed ratio 


9 

'“G 


5X, 


X, 


for X, 


X, 


9 

ax ^9X2 


X 


5/2 


for 


1X1 1 , 1X2 : 


and therefore inside the wake and at large distances or at a fixed X 


station and for large 'x^ 


5^G 

^XjClX^ 


is of smaller order than 


'ox 2 


1 


Since vj and are of the same order of magnitude, this implies 

that it is sufficient to examine only the rate of convergence of the 


integral containing 


o^G 


ax^2 


We denote by I that integral; 


I = 2p 


JJ 


bX. 


Y V^ v| dX^ cix., 
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whitli can also be written as: 


where: 


20d r +5d 


I - 2p 


r J r 




2d 


-5d 


2V,v'dX2 


!>dX^ + Il+l2 


+ “ r +5d 




h = ^v^vidX^^dX^ 


20d l-5d 


(3.23) 


(3.24) 


4- 00^ -f. 


I2 = 2p 


B^G 


Vi v[ dX^I 


2d I 5d ^^1 


+<»r-5d 2 

+ 2p J Jj ^^v'dx ldx , 
^ L 1-- ''31 


(3.25) 


and we propose to estimate the values of and . Introducing the 

non-dimensional variables as before and letting v^ , cp^ , cp^ , , Y.^ , 

, and denote these new variables we can write: 

.2 rfr / »2-^2>'-gr^i>' 1 1 - 




cp cos 2 tt( 0 -ft)+ cp cos 4TT(0^-ft) ) dx > dx 


or since cp^ is an odd function of X^ , cp^ '^i even functions 

of x^ : 


^1 . 2 r J r r ^^ 2 '^2 

..2 TtJ |J| 




(X,-YJ -(X^-Yj^) v“2'*2^ 


(X-+YJ^-(X,-Y, f 


"^2o|o \[(Xj-Yp2+(X2-Y2)^] [(X^ -Yp^+(X2+Y2)^3 


Z^\ ViTiC0s2TT(e^-ft)dX^dXj 


H-®r+ 5 /- 2 2 

^ (X^-Yp^-CXj-Y^)^ 



2 2 
^2 -^1 




[(X^-Y^)2+(X2-Y2)2] (x2+X^2) 


^ cp^ 4TT(92-ft)dX2 


dX, 
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5 2 _^ 2 




■ ^ J 1 J vj T2 4TT(9^_.ft)dX, I dX^ 


20 I 0 ^^1 "^^ 2 "^^ 


( ■^ , ?.h I 


As is clear from (3.26) and from the expression of in Lho form 

the shedding frequency component of the surface-pressure is the result: of 
an incomplete cancellation of the contributions coming from the parts of 
the wake below and above the X^-axis. This is best seen by putting ” 0 
in the above expression or in (3.19), the resulting contribution then 
vanishes. This is due to the fact that at a given time t an-l at fixed 
station the shedding frequency component of v| for instance, is 
antisymmetric with respect to the wake axis, and consequently the induced 
pressure at a point on that axis is the sum of two exactly opposite con- 
tributions; if, however, is not zero the cancellation is not complete* 

and there is a net contribution. The double frequency component on the 
other hand results from two additive contributions and is not the result 
of an incomplete cancellation effect. With this in mind it is clear that 
contributions from distant points in the wake are small not only because 
the velocity fluctuations at the shedding frequency there are small but 
also because the difference in distance between a point on the cylinder 
and two points sjrmmetric with respect to the X^-axis becomes small. This 
however does not apply when the double frequency component is involved. 

For these reasons, and in order to achieve better estimates, was 

(and will be) written in the above form. We now write: 




pu 


pu pU, 


2 


(3.27) 


pu 
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“^1 ^^2 “^3 

where -k , r , and k- are respectively the first, second, and 

pl'oo pUoo 

third term on the right-hand side of (3.26). An estimate of each of these 
terms as well as of \ ^ is found in Appendix B. The results are as 
follows: 


I 



0.0788X10"^ 


(3.28) 



^ 0.0190X10"^ 


(3.29) 



(3.30) 


(3.31) 


We note that the above estimates constitute only upper bounds to the 

Jl J2 ^2 

values of the integrals 5 * , r , ^ , and 5 - . These upper 

PUOD^ PUOD^ pUoo PUOD 

bounds are sometimes too gross (in view of the limitations of the formal 
procedure) to yield close estimates to the real values of the integrals. 

If we consider for instance the component of the pressure at the shedding 
frequency f , an estimate of the error involved in integrating I (see 
(3.23)) over the finite domain 2 ^ 20 , -5 - 5 is given by (3.28) 

and (3.31). By referring to Figure 1 we can see that an upper bound to 
that error is about 17 (given by (3.28)) plus 107 (given by (3.31)) of the 
maximum calculated value. The estimate given by (3.31) is too gross as 
is discussed in Appendix B and we can conclude, after assuming that a similar 
accuracy is obtained from the second term in (3.10), that to a very good 
approximation the contribution to the fluctuating surface-pressure at the 
shedding frequency comes from the near-wake region defined by 2 ^X ^^20 , 
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-5 ‘ 5 . On the other hand if we consider the component of Hu m \surt* 

at the double frequency 2f the estimates (3.29) tand (3.30) \d an upper 
bound for the error involved which is about 207 of the maximum calculatetl 
value. This, however, does not give any indication of the accuracv of the 
results which, we believe, is much better than is revealed by the above 
figure; and as for the shedding frequency component, we can conclude that 
the component of the surface-pressure at the dcvible frequency is determined 
to a large extent by the velocity fluctuations in the near-wake region. 

3 .3.^ Comparison with Other Results and Discussion . 

As we have said in Section 3.2, tlu only available experimental re- 
sults on the unst(*ady loading on the cylinder in the stable Reynolds num- 
ber range are those of Tanida et al. (1973). They measured directly the 
fluctuating lift force, but found that the fluctuating drag was too small 

I 

Cj j of the lift 

coefficient have been introduced already in Section 3.2. Comparison of 
our result (3.21) (at Re = 56 ) with these data shows excellent agreement 
and confirms the rather low level of lift fluctuation inherent to the low 
Reynolds number regime. By contrast the value for jc [ obtained by 

1 i 

Phillips (1956) at Re = 56 is in strong disagreement with the above 
results. This value ~ C.76) is about 40 times larger than ours 

(jc" 1 =0.0186) . Phillips used the same data of Kovasznay which have been 

Li 

used in the previous sections, and calculated the fluctuations in lift 
and drag according to the formulas; 
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pdT 




However, since behind an infinitely long cylinder the wake remains un- 
steady (if it is so near the cylinder) 
at arbitrarily large distances down- 
stream the above formulas are incom- 
plete and an additional non-vanishing 
term should be added to the right-hand 
vSide of each of the relations (3.31^) 
and (3*33). This can be seen by a]>- 
plying the momentum theorem to the 
fluid inside the control volume, hav- 
ing as cross -‘Unction and one dia- 

meter in depth in tlu* spanwise direction, 
in thi> fol lowing manner: 
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or after separation of the mean and fluctuating parts of each variable: 
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Taking the time average of the above relation and subtracting the result 


from the relation itself we obtain: 
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If now we designate by the ith component of the fluctuating force 

per unit span on the cylinder the above relation becomes: 
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On C. , R being very large, !vl'; er.cept on the v.’ake portion 
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where | vl 1 
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p' but since the flow is only npproximately t\;o- 


dimensional we cm assume fo"*- the present purpt^sc^ that it ]'• of smallt'r 
1 \ 

order than - . Since W~R‘ , upon increasing R indefinitely and iu;- 
suming that the first term on the right-hand side of (3,34) converge?- v.e 
obtain: 


F! - 
1. 


i '^(oV)dX,, dX - ^ p(v. vl 4* v.vl )n . ds 
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( 1 , 33 ) 


On W ^2 “ 1 ~ ’ assuming that v^ vanishes in the far-walu^ 

and taking i equal to 1 and 2 separately we obtain the follovnng re- 
la tions : 
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(3,3b) 
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(■3.37) 


^ - p ^ Jj J 

>: w 

(3.36) and (3.37) are to be compared with (3.32) and (3.33). In addition, 
we note that even in the above form it is doubtful whether (3.36) and 

(3.37) can be used efficiently to evaluate L and L' , for the converserce 
of the first term on the right-hand side of each of these relations is at 
best too slow and the integration would require an infinite amount of data, 
in fact, Phillips used the form (3.33) for the drag only, while for the 
lift he used the fact that: 

f'f’ f(* "^v' 

T. E ' ’ 

But since only v| is known from measurements he was led to differentiate 
the numerical data several times (for instance, cp^ was differentiated 
numerically three times with respect to ) to evaluate L , This, as 
we have mentioned earlier, can destroy the accuracy of the computations 
and lead to results totally different from the sought ones. 

Experimental results on the surface -pressure distribution when Re 
is in the stable range a'"e not available. Even at higher Reynolds numbers 
there are few instances where such results have been reported in the lit- 
erature. Those of Gerrard (1961) and McGregor (1957) are given in terms 

of intensity of fluctuations at the shedding frequency and its first har- 

4 

monic, at several values of Re between 10 and a value slightly above 

10^ . These two components (i.c., the shedding frequency component and 

its first harmonic) were found to be domrriant in the range studied 
3 5 

(4/10 Re 10 in Gerrard^s case). We reproduce below the results in 
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question as they were presonted in Figures 8 and 9 of Gerrard^s paper 
which also include the data of McGregor, The main feature of these re- 
sults is the invariance of the angular distribution with Reynolds number: 



REPRODUCIBILITY OP THi- 
ORIGINAL IS 


ict UK 8, Anpfiiliir distril>iitiou nf int^‘usi(y nf prf*sHun> ul tin* lundHtncntMl 

fmiiu*ri(’y. loUioE: f i. 41; A, 4-;i; X , o, r>-2; 4 , McHiv^nr 4*04* 



Anp:ular jmsilirai 0 H 

I'lt.niK 0. Angular din(n'bulion (»f int<*Msily (4’ surftu'c’ jirr^sura at tlio .second harrnnnit^ 

I'h'fjuoncy. 1 jn. diaTii. f\\liiid(T; x , 3 in. diam. fylinder; 4 , l^IcCrrgnr. liin. fliain. 

the intensity of fluctuations at the shedding frequency is always r.aximum 

at about 120^ from the front stagnation point and drops substantially at 

the front and back of the cylinder, while that at the double frequency is 

maximum at the back of the cylinder. In addition, Gerrard observed that 
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in the range studied the pressure is essentially in phase over one side 
(upper or lower) of the cylinder and 180° out of phase with that on tlie 
other side. Referring back to Figures 1 and 2 we can see that a remark- 
able similarity exists between our results at low Reynolds number and the 
above ones at much higher Reynolds numbers. The question whether the in- 
variance of the angular distribution with Re can be extended to the 
stable range may tentatively be asked. In order to facilitate the com- 
parison, the results of Figure 1 were plotted again, in the manner fol- 
lowed by Gerrard in Figures 8 and 9, and are shown in Figures 3 and 4. It 
readily observed that the two sets of figures are similar to a great de- 
free, the difference being that at high Reynolds numbers, the intensity 
of fluctuations at the shedding frequency drops more slowly when the 
front and back stagnation points are approached and the ratio of the in- 
tensity of fluctuations at the double frequency on the back of the cy- 
linder to that on the front is greater than it is at low Re . We defer 
the discussion of these similarities and differences to the next section 
where a general discussion of the high Reynolds number case is given. 

To conclude the present section we examine the effect of varying Re 
in the range below Re = 90 . As we have said at the beginning of 3,3, 
the nature of the velocity fluctuations in the wake does not change when 
Re is varied from 40 to some value below 90 : the fluctuations de- 
velop always as a result of the instability of the laminar wake and the 
flow remains two-dimensional. The only change is in the intensity of 
fluctuations and in the distance from the cylinder at which this intensity 
reaches its peak value. Upon increasing Re the velocity profiles in the 
laminar wake become steeper, the concentration of vorticity in the street 

becomes stronger, and the distance at which the intensity of fluctuations 
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Figure 3. 
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reaches a maximum value becomes shorter. Consequently, the pressure fluc- 
tuations on the cylinder, at the shedding frequency, for instance, continiu 
to be the result of an incomplete cancellation of two contributions coming 
from the parts of the wake below and above the centerline (see ; 3.3.4) 
and the magnitude of the fluctuating lift can be expected to remain small 
even though it tends to increase slightly with Re as a result of the in- 
crease in the intensity of the velocity fluctuations in the proximity of 
the cylinder. Similar increase may be expected for the fluctuation in 
drag for the same reasons. 

3.4 General Discussion - Reynolds Number Effect . 

When Re exceeds some value near 90 , the flow becomes sensitive 
to small three-dimensional disturbances or slight non-uniformities in the 
cylinder or in the flow upstream, and some of the details of the wake pat- 
tern become uncertain and dependent upon the nature of these disturbances 
and non -uniform! ties, in addition to their dependence on the details of 
the experimental setup, as has been already discussed in Section 3.2. 

With this in mind it is clear that it would be impossible to attach a 
unique value for the magnitude of the fluctuating lift or drag coefficient 
to each value of the Reynolds number Re . This is true whether we are 
talking about laboratory experiments or about real life situations as they 
occur in applications. The theory of the preceeding chapter, on the other 
hand, allows certain sys temization in the study of the fluctuating loads, 
and this by identifying the role of the cylinder geometry and that of the 
various components of the wake flow, as they may occur, in bringing about 
such loads. On the basis of that theory and o ir experience with the low 
Reynolds number range treated in the preceding section, we will attempt 
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bolow to speculate on, or when it is possible, "predict" the maimer in 
which the unsteady loading on the cylinder would vary with Re under var- 
ious conditions. At the same time some comments on previous attempts will 
be made and areas where further research is needed will be indicated. 

Using (2.39) and (2.40) the general solution (2.15) given in the pre- 
ceding chanter can be written in the form: 


P(Y,k,t) 



a^G(xjY;k) 

c^X.^X, 

1 J 


{2v,V, +F[v!v'. 
1 .1 1 J 
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-f 2ipk JjG(XlY;k)C(v^:2^ ' 

r. 

- pk“ JjG(XlY;k)[2(v^V^ +V2V2)+F[v|^ + vV']]dX 
E 

-t M ^G(X1 Y;k)n.V“Vds . (3.38) 

I 

C 

From the above relation it is clear that, insofar as the pressure P is 
concerned the two aspects of the fluctuating wake, namely the two-dimen- 
sional one represented in (3.38) by the first, third, and fourth integral 
and the three-dimensional one represented by the second Integral can be 
examined separately, and both play an important role in the determination 
of P . Some important features of these two aspects are examined below 
in detail . 

Considering first the two-dimensional aspect, i.e., the wake flow as 
it appears in the plane (X^ , X^) and concentrating on the vortex street 
alone, one can argue that the intensity of the velocity fluctuations Vj, 
and and the manner in which it (the intensity) changes with Re 
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should have an important role in determining the magnitude of P ami its 
variation with Reynolds number. The intensity of the velocity fluctuations 
on the other hand is closely related to the strength of the vortices in 
the street. Measuring the strength of these vortices is a difficult pro- 
cess which is usually accomplished by matching a theoretical model of the 
flow to hot-wire measurements in the wake; this gives only an average 
value and has been attempted by several Investigators. It has always been 
considered as remarkable that only a part of the rate of vorticity gener- 
ated in the boundary layers on the cylinder is found in the vortex street 
(see Berger and Wille 1972). If we define (as in Berger and Wille\s re- 
view) the rate of discharge of vorticity in one row of the street by: 

K = r f = Au ^ , 

'JO ’ 

where F is the circulation of a vortex, and the initial rate of discharge 
of vorticity from the boundary layer on one side of the cylinder by: 

K = 0.5U " 
s s ’ 

where is the velocity at the outer side of the boundary layer near 

the separation point, the ratio: 

K 2A 



where C is the base-pressure coefficient is a measure of the fraction 
Ps 

of the vorticity which escapes annihilation during the interaction of the 
shear layers from both sides of the cylinder. This ratio is found to be 
about 0,4 in the low Reynolds number range corresponding to Tritton's low- 
speed mode and to vary between 0.5 and 0,6 or 0.66 at higher Reynolds 


87 


numbers (see Berger and Wille 1972). On the other hand a measure ol‘ tlii' 
strength of the individual vortices in the street is F . From the above 
relations it follows that: 

1-C 

J1 ^ A !k A 

U d St 2St K 

00 5 

C varies with Re ; it is about -0,5 when the wake oscillations begin, 
^ s 

it decreases to about -0.9 at Re = 200 , then increases to about -0.7 

3 5 

at Re = 10 , and finally decreases to -1.2 near Re - 10 before in- 

creasing again at the critical Reynolds number (see Roshko and Fiszdon 
1967). As for St , it is about 0.12 at the onset of oscillations, it 
increases to about 0.19 at Re = 200 , then to about 0.21 at Re = 10^ 
and remains near this value for all subcritical Reynolds numbers (see 
Roshko 1954). After some calculations one can easily see that the ratio 

(1-C )/2St varies little in the range considered; it takes the values 6, 

"" 3 S 

5, 4, and 5 at Re = 50, 200, 10 , and 10 respectively. In view of 

the small changes of with Re we note here that it is even more 

remarkable that the strength of the vortices in the street (non-dimension- 
alized) varies little over the whole range of Re below the critical 
value. With this in mind, if we look at the results of Tanida et al. 

(1973) and those shown in Figure 15 of Morkovin (1964) (both introduced 
in Section 3.2) concerning the magnitude of the fluctuating lift coeffi- 
cient we see that this coefficient increases from very low values in the 
range below Re = 100 to high values near Re = 10^ and this despite the 
increasing three-dimensionality of the vortex filaments in the street with 
Re . Therefore, neglecting for the moment whatever contribution to P 
the flow-wise vorticity term (the second on the right-hand side of (3,38)) 
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has, the absolute strength of the individual vortices In the street and 
the closely related intensity of velocity fluctuations there are not of 
prime importance in determining the magnitude, of P , What is important 
rather is the net strength of the total vorticity (mean and fluctuating 
part) in the wake, which does not vanish at all times as in the case of 
an idealized Karman street, and upon which the intensity of velocity flue 
tuations in the vicinity of the cylinder depends to a large extent. This 
can be seen In the following way. Let us consider a cylinder started from 
rest in an infinite fluid equally at rest. As a result of the no-slip 
condition, vorticity is generated at the coundary and is mainly swept to 
tiie back of the cylinder (at sufficient high Re ) to form the wake. By 
applying Kelvin's theorem as in airfoil theory it is seen that an equal 
amount of vorticity must be placed inside the cylinder so that the circu- 
lation on a contour surrounding the cylinder and the part of the fluid 
which has already recei ;ed some vorticity is zero. Before the onset of 
instability the vorticity generated at the boundary is antis 3 mimetrically 
distributed with respect to the wake centerline, therefore the net vorti- 
city in the wake vanishes and so does the bound vorticity inside the cyl- 
inder, This implies that no circulation exists on a circuit closely sur- 
rounding the cylinder and consequently the cylinder does not experience 
any lift force. On the other hand, when high shear region in the wake 
becomes unstable to small disturbances, antisymmetric (with respect to 
the wake centerline) n the longitudinal velocity, modes of disturbance 
vorticity, symmetric with respect to the wake centerline, grov? and travel 
downstream in a wavelike pattern (see Mattingly and Crhninalr-. 1972). This 

An idealized Karman street, if left alone in an infinite fluid, does not 
move in a transversal direction. 
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additional vorticity is simply sui>erposed on ihu original mean veloeity 


distribution in the early stage of its growth while further dowustrevim it 
interacts with the mean vorticity, modifies it, gets modified itself and 
Degins to diffuse In the fluid. However, at all times the mean vorticity 
remains ant iysmme trie with respect t(^ the wake centerline and the distur- 
bance vorticity at the fundamental frequency remains symmetric. Therefore 
the net vorticity in the wake diH*s not vanish and alternates (at the fun- 
damental frequency) between positive and negative values. In addition, by 
the nature of the process of generation 
of the disturbance vorticity, tiils net 
strength is small; if for instance at 
a given time, and at each downstream 
station v;e sum the vorticity in the 
transversal direction and plot the 
obtained values as a function of the downstream coordinate we obtain 

a curve of the form shoxm in the figure, the net strength in question is 
nothing but the algebraic area under this curve, this is of the order of 
the area under the curve AB and therefore is small. This implies that 
the bound vorcicity is small and so is the lift force. As for the velocity 
fluctuations in the wake, they are due to the field of the unsteady vorti- 
city both in the wake and in the cylinder. Near the cylinder the field of 
the bound vorticity is dominant, and since in this case the latter is 
small, the velocity fluctuations there are small. 

What we have described above is the state of affairs corresponding to 
the low-speed mode of vortex shedding. This is to be contrasted to the 
case when the vortex shedding results from the direct interaction of the 



90 


two shonr hiyors emanating from the surface of the cylinder (correpomling 


to the high-speed mode of Trltton) the 
(19f56b) and observed by Mattingly (196 
stuay on the non-linear interaction of 
nathy and Kronauer (l'»62) applies. 
Without entering into details (re- 
ference is made to the mentioned 
papers) and referring to the figure, 
the net vorticity in the wake (as- 


way it was described by tlerrard 
1 ), and to which the fundamental 
two infinite vorcex sheets of Aber 



Sliming that all the vorticity is 
concentrated in the. street) is 

dominated by the difference in strength between vortex B and vortex A . 
At certain moments which occur periodically this difference is largo; the 
induced bound vorticity is then large and the lift has a maximum amplitude. 
In addition, the velocity fluctuations in the vicinity of the cylinder, 
which are due mainly .o the field of vortex B and the bound vortex, reach 
high values and are in phase in a direction parallel to the X^-axis. 

The three-dimensional nature of the flow in the wake, when Re is 
above 90 , is far less clear and much more difficult to describe. The 
pieces of information available are scattered in the literature and are 
often reported in a rather crude manner. For instance, spanwise corre- 
lation lengths determined using two hot-wires separated in the x^-direction. 


"This is the reason that calculated convection speeds, of the street vor- 
tices, based on phase measurements increase without limit in the vicinity 
of the cylinder (see Bloor and Gerrard 1966, and Simmon^ 1974). Simmons 
recognized the existence of the bound vortex, but assumed that the vor- 
tices in the wake are equivalent to a distribution of vortex doublets, 
while in fact it is the total vorticity including the bound vortex which 
is so. 


91 


are usiiallv given without spoci lying whether the signals from the hot- 
wires were filtered at the shedding freqmmcy or not. Similarly it j:- 
wot clear sometimes whether the perit^dic i ty in the spanwise direction, 
observed at certain values of Re , is associated with slantwise shedding 
or simply with waviness in the vortex filaments. Hven the more special- 
ised studies of clerrard (l*ih6a) and Gaster (1969, 1971) on the effect of 
spanwise iion-mu formi t les (in the cylinder or in the free-stream) leave 
some important issues unanswered. In addition, Cierrard (19o6a) attributed 
the existence of low f retjuency modu I a i Ions , in the transition range, to 
the slmul tansous existeiu e in a direct ion along the span of laminar and 
tiirbulvnit vortices, which also resulted in the continuation of the vor- 
tex filam('nts into the direction of the free -stream. He also noted that 
in the irregular range, when all the vortices are turbulent, this effect 
is less pronounced. However, he later showed (Gerrard 1967b) that suppres- 
sing the vortex strc'et with a splitter plate did not alter the low fre- 
quency coupHuient in the vicinity of the cylinder, thus proving that this 
component is unrelated to the vortices in the street. 

Leaving aside such factors as the spanwise non-uniformity and the details 
of the experimental setup, one thing we can be sure about is that the pre- 
sence of three-dimensional disturbances in the free-stream will ultimately 
render the flow in the near wake three-dimensional. In Section 2.6 we 
have suggested that such a phenomenon takes place in the individual shear 
layers, and gave as evidence the observations of Humphrey (1960) at the 
critical Reynolds number. Other pieces of evidence can be given to support 
further this point: for instance, the observation of Hama (1957) that 

’’the ransverse waviness of the vortex lines appears almost immediately, 
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In addition, wt‘ cnn nr- 


turtlier develops apparontly in the shear laye»r." 
gne that a f roo shear layer is more likely to be vulnerable to small dis- 
turbances than a rolling or rolled up vortex. It is also tempting to 
conclude, ns we did in Section 2.6, that as n result of the growth of 
these disturbances, flow-wise vortieity is generated, which either takes 
the form ol wnviness in the vortex filaments, or causes the slantwise 
shedding of vortices, or even leads to the low frequency Irregularities. 

If we accept this sequence of cause and effect, the various observations 
mentioned above become less puzzling. Another item of observation seems 
to fit the above picture: by oscillating the cylinder at the shedding 

frequency, Berger Tsee Wille 1966) has shown that laminar vortex streets 
with naralh^l tas opposed to slantx^;ise) shedding could be upheld until 
Re 350 . It is possible that upon oscillating the cylinder the separated 
shear layers become more stable against both small and large scale three- 
dimensional distruhances , thereby delaying the transition to turbulence 
and preventing the formation of flow-wise vortieity which v;ould otherwisr 
lead to slantwise shedding. Let us emphasize that: the preceding dlscussliui 
is merely speculative and is not intended to be definitive. Nevertheless, 
as is clear from the general solution (3. '38), the knowledge of the three- 
dimensional nature of the flow in the near-wake or the form under which it 
manifests itself under various conditions is essential for a successful pre^- 
diction of the unsteady loading on the cylinder. This, we believe, can bo 
achieved only by conducting experiments under controlled conditions. 

We now proceed to the discussion of the changes in the magnitude of 

Slantwise shedding is equivalent, from a kincmatical point of view, to a 
vortical component in the direction of the mean flow plus a certain perio- 
dicity in the spanwise direction. 
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the unsteady loading with Reynolds numln*r. The lirst question one is 
faced with, when Re exceeds 9(J , is which of the two modes of vortex 
shedding discussed above is likely to occur. Is it the low-speed mode 
or is it the high-speed mode, v^e know that the low-speed mode is caused 
by the instability of the laminar wake and is triggered by the presence 
of lineari^ably small di st urbanct'^: . On the other hand, the high-speed 
mode can imly be trigge*red by I’initi' disturbances (Abernathy and Krou- 
aner Lf we increase Re ^*,iadually from values below 90 wherc^ 

the low-speed modi* is prevailing, tlu? velocity fluctuations induced in the 
vicinity of the cylinder, which are antisymmeLric in the longitudinal com- 
ponent and thereiore constitute a potentially triggering agent for the 
high-s]H*ed mode, incrtsn-U‘ but remain very small at low Re and therefore 
are not able to caiuu* the iioxi-linear interaction of the two separted shear 
layers. If on the otlier hand the level of disturbances in the free-stream 
is high enough, and Re is increased beyond 90 the absence cf the two 
standing eddies (l.e., the fact that the two shear layers spring freely 
from the sides of the cylinder) will lead to this type of interaction. We 
don't know the background turbulence level in Tritton's experiment of 1959, 
however in the 1971 paper Tritton mentions that this level was high, Ko- 
vasznay (1949) who had a background turbulence level of 0.06X did not 
observe the transition of Tritton and in his words, "the vortices develop 
some distance downstream within the Reynolds number range 40 to 160 , 
and are not shed direct'^y from the cylinder, conseque:itly the phenomenon 
can properly be considered as an instability of the laminar v;ake, that 
develops to an amplitude limit but dies out before becoming turbulent," 
Therefore, we can conclude that with low free -stream turbulence level the 
wake remains in the low-speed mode even after Re = 90 is exceeded; the 
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increase in the lift and drag fluctuations with Re is then moderate and 
is quickly counteracted by the increasing effect of three-dimensionality 
with Re (a typical value of the spanwise correlation length for 
90<Re< 150 - 200 is 17d , see Phillips 1956). Referring back to the 
dcj ca of lanida et al. (1973) in the stable range, which were taken in an 
oil tank where presumably che level of disturbances is very low, it is 
very likely that the measured values of 1 correspond to the low-speed 

ij 

mode of vortex shedding. On the other hand the relatively high value 
(ic*^! =0.27) obtained by Jordan and Fromm (1972) at Re = 100 might have 

Lj 

resulted from prematurely triggering the high-speed mode of vortex shed- 
ding, by twisting the cylinder back and forth In order to perturb the 
numerical solution. 

When Re is increased beyond the stable range, the increasing effect 
of three-dimensionality corresponding to the second term in (3.38) becomes 
more difficult to assess (a typical value of the spanwise correlation 
length in the transition range is lOd given by Eoshko (1954) at 
Re “• 220) . However, we expect that, like the Strouhal number, the lift 
and drag fluctuations are not well-defined in the transition range, in 
particular when Re is close to 400 (see Section 3.2). For this and 
other reasons it is doubtful that the value jc j =0.75 obtained by 

Ij 

Jordan and Fromm at Re = 400 can be considered as typical of the range 
in question. 


For a given two-dimensional distribution of velocity fluctuations in the 
wake, the difference in magnitude of he pressure fluctuations corre- 
sponding to an infinite or a finite spanwise length scale 2rT/k is the 
same as the difference in temperature distribution in a medium corre- ^ 
spending respectively to a zero or a finite coefficient of absorption k 
(see Section 2.2) • 


95 


The possibility of a coupling between the transition waves of Hloor 
(1964) and the shedding frequency mode extends to Re = 1.3X10^ . To our 
knowledge no measurements of unsteady forces have ever been reported in 
this range, and one may well wonder whether this is due to measurement 
difficulties of the type encountered in the transition range. Above 

3 

Re =1.3x10 the transition waves and the shedding frequency component are 

decoupled, and therefore the lift and drag fluctuations at the sheeding 

frequency and its first harmonic respectively can be expected to be well 

defined as in the stable range. However, the qu stion concerning the mode 

of vortex shedding arises aga in. In addition, the wake becomes highly 

three-dimensional; typical values of the spanwise correlation length are 

5 

three to five or six diameters up to Re = 10 , and we expect the contri- 

bution from the flow -wise vortical component to become important. Conse- 
quently the changes in the magriitude of the nmstcSdy loading with Re in 
this range are the result of the simultaneous change in both the two-di- 
mensional and the three-dimensional characteristics of the flow in the 
near-wake. Since very little is known about the three-dimensional char- 
acteristics, it is very difficult at this stage to speculate on their 
possible effect. On the other hand the two-dimensional characteristics 
and their effect on the magnitude of the fluctuating lift were the subject 
of a series of three papers by Gerrard (1965, 1966b, 1967a). Since then 
two relevant papers appeared in the literature and we feel that some com- 
ments here are in order to further emphasize some uncertain issues. 

On the basis of a potential flow model Gerrard (1967a) suggested the 
possibility that, when the turbulence level in the free-stream is very low, 
the wake remains in the low-speed mode up to values of Re of the order 
of a few thousands. However he later questioned that possibility, because 
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of the three-dimensional nature of the flow at such high Kcynolds nimibi r- 
and the invariance of the Strouhal number, which normally can be expected 
(the Strouhal number) to depend on the mode of shedding. In their photo- 
graphic study, made in a mercury tank, Papailllou and Kykoudis (19741 

3 

observed a wake in the low-speed mode at Re = 3.4>10 and a formation 

region whose length is in agreement with the value measured by Bloor 

(1964). This could have provided a confirmation of Gerrard's suggestion 

if it was not for the fact that the only picture showing a symmetrical 

formation region was taken with a cylinder v'hose span to diameter ratio 

is only 2.5 . This we know can impose the two -dimensionality on the flow, 

at least in the near-wake region. In addition, the recent data of Tanida 

3 4 

et al (1973) in the range 2X10 <Re<10 , which (the data) correspond 

to a low free-stream turbulence level, do not agree with those of Gerrard. 
They do not agree with those of TCeefe (1961) either, but Keefe had a 
background turbulence level of 0,3Z . On the other hand, the similarity 
of the uata of Gerrard and those of Keefe with the clearance holes open, 
suggests that perhaps the results of Gerrard suffer from end effects. If 
this is true, then the reason for the unrealistically large formation re- 
gion generated by Gerrard (1967a) numerical model, when forced to yield 
the very small fluctuating lift he measured, becomes clear, Gerrard 
(1967a) argued that with a hot-wire in the flow it may be impossible to 
obtain a S 3 mimetrical formation region, hence the discrepancy with Bloor 
(1964) is natural. But Papailliou and Lykoudis did not have a hot-x^ire in 
the flow and yet the length of the symmetrical formation region they ob- 
served is in agreementwith Bloor *s measured value. Another factor, which 
might have contributed to the large formation region in Gerrard *s model, 
is the fact that this model is purely two-dimensional. 
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To conclude, we consider again the question (noted in Seition 3.3,5) 
of the invariance with Reynolds number of the angular distributior of the 
surface-pressure fluctuations at the shedding frequency and its first 
harmonic. The similarities in this respect, between the stable range and 
the subcritical range, can be attributed to the fact that whenever a de- 
finite vortex shedding occurs at high Reynolds number, as is the case for 
instance in the subcritical range, the two-dimensional characteristics of 
the velocity fluctuations at the shedding frequency and its first harmonic 
are similar to those in the stable range. For instance, the components 
Vjj^ and v ^2 are antisjmimetric with respect to the wake centerline, and 
the components Vj 2 ''^21 symmetric. From this and the continuity 

equation one can easily deduce that and (the second subscript 

being used as for v| and v^ ) are symmetric with respect to the wake 

mm*- ' * 

m * 

centerline, and co ^2 ^21 antisymmetric. Referring to (3.38) 

and noting that G(XlY;k) is simply a geometrical function which depends 
only on the spanwise length scale and, like G(x1y) , on !x-y 1 , we can 
see upon inspecting the different integrals (concentrating only on the 
linear terms) that, as in the .:wo -dimensional case, the component of P 
at the shedding frequency is the result of an incomplete cancellation of 
two contributions coming from the regions below and above the wake center - 
line, while the component at the double frequency is the sum of two such 
contributions. This implies that the angular distributions of the surface - 
pressure at these two frequencies should have roughly the same form as at 
low Reynolds number. The difference between the two cases, on the other 
hand, lies in the fact that, when is put equal to zero (corresponding 

to the front and back stagnation points), the cancellation is not complete 
for the component of P at the shedding frequency, as it is at low Re . 
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This is due to the fact that at high Reynolds numbers the wake ’physical' 
centerline does not coincide at all times with the wake 'geometrical' 
centerline (i.e., the Xj^-axis), This phenomenon is often referred to as 
the wobbling of the wake axis and is manifested by a non -vanishing signal , 
at the shedding frequency, from a hot-wire placed on the X^^-axis isee 
Bloor and Gerrard 1966), As for the double frequency component, the fact 
that k is finite at high Reynolds numbers implies that the velocity 
fluctuations at the double frequency in the immediate vicinity of the cyl- 
inder supply a more dominant contribution to the surface -pressure than they 
do at low Re . This follows from the asymptotic behavior of G at large 
values of the argument (see Chapter 2) and the effect is more pronounced 
the larger k is. Hence, the ratio of the distance between the back 
stagnaLior, point and the rtigion supplying the dominant contribution at 
that point to the distance between the front stagnation point and the cor- 
responding region which supplies the aominant contribution there is smaller 
at high Reynolds number than it is at low Reynolds number. This may be 
the source of dissimilarity between our Figure 4 (at low Re ) and Gerrard 's 
(1961) Figure 9 (at high Re ) . 
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APPENDIX A 


The Role of the Nonlinear Terms In the Source Distribution Function 

The source distribution function in its original form, i.e,, as given 
by (2.5) and (2.6), can be written as: 

r^v 5v' Bvl <)v'. 9v' Bv' 

• 2p p - p 

1 i .1 i 

The first term on the rig’^ t-hand j*ide of (A.l) represents the amplifying 
effect of the mean shear on the velocity fluctuations, while the remaining 
two represent the effect of these fluctuations on themselves. In all the 
previous studies which considered the pressure fluctuations on a plane sur- 
face supporting a turbulent boundary layer (Kraichnan 1956b, Lilley and 
Hodgson 1960, and Lilley 1963), the latter effect was neglected. The rea- 
son being that duo to the large mean shear in the boundary layer, the con- 
tribution from the quadratic terms is only a fraction (4 to 67.) of the 
total contribution. However Corcos (1964) argued that the linear term 
(the first on the right-hand side of (A.l)) which is large only in a region 
very near the wall supplies a dominant contribution only at high frequen- 
cies, while at low frequencies the remaining terms have contributions at 
least comparable with that of the former and their scales and convective 
speed are more typical of those of the observed wall pressure. Nonethe- 
less Panton and Lineberger (1974) in their recent study retained only the 
turbulence -mean shear teirm and found that their results were in good agree- 
ment with observations for wavenumbers (in the streamwise direction) 

as low as , where 6 is the boundary layer thickness. 

In the wake behind a bluff cylinder the situation is different. The 
shear layers which spring freely from the sides of the cylinder increase 
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in width with the downstream dis- 


tance and at a faster rate the 
higher the Reynoldj^ number is. 

Although the velocity fluctua- 
tions do not peak very near 
the boundary where the mean shear 
is highest, they do peak, at each 
downstream station , at the 
center of the shear layers. This 
is true not only for the high fre 
quency components which result 
from instability of the individual shear layers themselves but also true 
for the shedding freqjuency comp :ment^ (see GerrarjJ 1967b). Across each*- 
layer the change in the mean velocity is of the order of the free -stream 
velocity while the change in the level of velocity fluctuations is 

about one order of magnitude smaller (see Figs, 2 and 3 of Hanson and 
Richardson 1968, where contours of constant mean velocity and r.m.s, 
values of velocity fluctuations in the near-wake are shown at two differ- 
ent values of Re : 10600 and 53000). With this in mind, the first term 

on the right-hand side of (A.l) can be expected to yield a contribution 
about one order of magnitude larger than that of the remaining tern.s. 

Note however that at values of the wavenumber k (in the spanwise direc- 
tion) or the frequency co (in time) for which the velocity fluctuations 
in the wake are not particularly significant. The linear term in q is 
not dominant and the main contribution to the fluctuating pressure at these 
wavenumbers and frequencies comes from the quadratic terms. 
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For the particular case of Reynolds number equal to 56 which is 
treated in detail in Chapter 3, the same argument given above is still 
valid. Inspection of the mean velocity profiles and those of the velocity 
fluctuations shows that at each downstream station the peak of the fluctu- 
ations occurs in the high shear region and that the gradient of these 
fluctuations in the transversal direction is at least one order of magni- 
tude smaller than that of the mean velocity. Therefore the contribution 
from the nonlinear terms in q can be expected not to exceed 107, of the 
total contribution. On the other hand, and as was noted above, the compo- 
nents of the pressure fluctuations at the frequencies 3f and 4f can be 
accounted for only by considering the nonlinear terms. But these compo- 
nents are negligibly small and practically non-measurable at this low 
• Reynoljds number. • • • 
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APPENDIX B 

Estimation of Some Integrals 


Estimation of J-/dD . We first consider J, . We have; 
1 «> 1 


pu 


+ ”('+5 2 ^ 2 '^2 0 

20 Co ^ ' 


-t' 1 

l-V,, ] I X 


cos 2 TT( 0 j-ft)dX 2 > dx^ 


= .ill f jl- 

ttJ 1 .1 2 „ 9,„2 


20 V. 0 ^^2 


(X^-Y^)-+X 2 


-2X,Y,,+Y„‘' I 
] -iog[ 1+ — 2 ] f y 
(x^-Yi)-+X 2 j 


cp^ cos 2Tr(0j^-ft)dX^ f 

- 

1 1 

Since -“-Y^ have in the domain of integration: 


+ 2X2Y2 + Y2^ « (X^-Y^)Vx2^ 


By expanding the logarithms and retaining only the first term in each 
expansion, we obtain: 


T ■'■“r ^ 2 

•^1 ^ _i r Jr d 

TT 2 ttJ 1 tl 2 

20 I 0 ^^2 


20 ^0 5 X 2 (Xi-Y2)+X2 


‘V 2 — 

r^3vj <Pj cos 2 n(Sj.£t)dX 2 I dXj 


^^2 r J r C 

“J 1 J ~T 


[ 


-X 2 


20 V. 0 s^2 +^2 


^1 2n(9^-ft)dX2 j dX^ 


or, 


^1 "^^2 r"” 

— ^ = ~ I K^(X^) cos 2TT(9^-ft)dX^ 
20 


(B.n 
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where: 


K, (X/) = 

11 K 


X,“ (Xj-Y^',"+x,; 


V, dx, , 


— "^2 

in which v, , , and -.7 [ 7^ arc all positive functions 

(X -Y )*-+x,,^ 

of X,, for Xj '* 20 . Therefore, k|(X|)‘ is a positive function of Xj 

for X > 20 , In addition, K^(X^) is a monot lically decreasing 

function of X for X, '*20 , because 'n. and 7; [ ^7 7;] are 

_ ^ ‘ ^ =»X2“ (Xj-Yi)-+X2“ 

while Vj is a slightly increasing function of X^ . An upper boimd 

can be four’ for Kj in the Intel il 20 - X^ lOO . By examining the 

data on the velocity fluctuations w can see that: 


• — , for 20 • X loo ^ 

1 


and since v^ 1 , we have (for !’o loo ) 


■1 f - -X') _ 

K.(X ) f- ' 7 [ % 

^1 P^X^^ (X^-Y^r+X^^ 


, f 25-(X -Y ) ^ ^ 

K (X ) r f-) -V--~2 + —7 ^ 

^1 |C(X -Y )2+;;5] (Xj^-Y )“^ 


2 3X25(X^-Y^)‘+(25)^ 

^1 [(Xj-Yp2+25]^(X^-Y^' 


2 3X25(X^-Y^)^+3X(25)^ 

^1 C(X^-Y^)^+25]^X^-Yp^ 
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0 


3X25 


^ 

[(X^-Yj)^+25](X^-Y^)^ 


<- A 3X25 
■ (X^-Y^)" 


and therefore: 


K^(X^) •: 


150 

Xj(Xj-Yj)' 


for 20 Xj •?• 100 . 


In estimating we will neglect the part of the integral beyond 

Xj = 100 assuming that It is negligibly small and does not affect the. 
order of ‘magn' tude of the retained jjart.' VJith this in mind we write 
instead of 05.1) 


''+Y,, 




100 
_2 f 

TT li 
20 


K^(X^) cos 2TT(6^-ft)dX^ 


(B.2) 


Since K^(X^) is a positive monotonically decreasing function, it follows 
from the second mean value theorem (see Whittaker and Watson 1963) that: 


J, 4Y^ f‘l 

— t- = — K, (20) cos 2tt(9, -ft)dX, 

,, TT 1 •' 1 i 

20 

where is such that 20 100 . Now we know that given a linear 

function «(X^) such that: 


In the far-wake cp^ decays more slowly than 1/X. but it is *v.asonable 
to assume that the velocity fluctuations beyond = 100 do not affect 
the surface-pressure field. 
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the following Integral: 
b 

f cos 2nQ! dX 

1 

a 


where n and b are arbitrary but finite, is bounded with -^5 as 

n/i 

upv>er bound, is an almost linear function of , for the present 

purpose it can be assumed linear and inspection of the data shows that: 



dX^ 15 


Therefore we have: 



-K,(20) 


TT X 


15 


or; 

^1 j _1_ 15 X 150 

pU^^! n“ 20(20-Yj)^ 

,, J_ 15 y 150 

‘2 4 ’ 

rr^ 20(19.5) 

and finally: 


1 

5" ■: 0.0788 X 10 

DU I 
00 


Estimation of 




We have: 


(B.3) 
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2 2 

Xi +X2" 




Again, since in the domain of integration: 


y/+y/- 2 XjYj- 2 X 2 Y 2 « X/+X/ 


by expanding the logarithm and reta-* ■- :.ag only the part of the first term 
which is of lowest order: 


J2 2 ^ Y,X,+Y„X 


PU„^ "" 20 I -5 X 


X +Y„X '] __ 

2 1 ''1 fi • 


or since is even in X2 , ^2 even in X2 , and since rp2 (by 

inspection the data) is different from zero only for - 3 :SX 2^+3 and 
X^ 100 : 


'^2 _ 
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ff 
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can be vjritten as: 
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J, _ 2Y^ , 


UK) 


pu 


20 


K2(X^)cos4TT(02-ft)dXj , 
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where: 


K2(X^) 


+3 , 

1 ^ . 

[ ■"'0 

-3 ^'' 2 ^ 

[xj^+x,"] 


Vi cp^dX^ 


As before, it can be shown that K^(X^) is a positive monotonically de- 
creasing function of . In addition, since (by inspection of the data) 

^ 20 - Xj 100 
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Kj(Xi) S 
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Now from (B.4) and in view of the second mean value theorem it follows 
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Kstimation of J<,/pU 
3 °° 

The cancellation effect mentioned in Section 3.3.4 cannot be ex- 
ploited to find an upper bound for . However, in order to improve 

the estimation we can use the fact that: 

2 

.. 0.7 ,, ^2 V 
'^2 ’ Xj 9 

We write: 
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is a i'ositive monotonical ly decreasing function of and: 


hi -h ix{-+x/J 


x^ 

(l--|-)dX2 
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1.4 , 3 ^ 3 

‘ T, 
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Also, we have as before: 
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0.1329 X 10 
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2 

Estimation of I /pU 
/ 

2 

Iz/pU^ accounts for the effect of neglecting the contribution from 
the ve^ icity fluctuations outside the strip jx^l - 5 . There, only the 
fluctuations at the shedding frequency are present and they are very small* 
Upon examining the data we find that: 
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„ .. 0.0175 .. „ c 

ffl, 0 X- ^ 5 , 


where again tne variables are non-dimensional. From (3.25) it follows 
tha t : 
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We note that: 


2 r ] 2X2[3(Xj-Yj)2.x^2^ 


2.„ 2 


^X, l(X^-Yp"+X 


> = 


[(Xj-Yp2+x^2]^ 
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and therefore in the domain of Integration (i.e., 2' , 

5 X.^ +*» ) this quantity changes its sign. It follows that the proce- 

dure used in the previous estimations cannot be applied here to account 
for the behavior of the cosine function in (B.7). Nevertheless, we will 
carry out the estimation below by taking the absolute values of the dif- 
fprr*xit terms in (B*7), knowing that only a gross estimate can emerge: 



Now, we have: 
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The maximum value of the term containing the logarithm on the right-hand 
side of the above inequality cannot be determined accurately, we therefore 
evaluate that term at ^ corresponding to the maximum value of 

(see Figure 1). The resulting estimate is: 



^ 0.6308 X 10 


(B.8) 
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